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ABSTRACT OF THE DISSERTATION 


Nonlinear Quantum Search 


by 


Thomas Giechaung Wong 
Doctor of Philosophy in Physics 
University of California, San Diego, 2014 
Professor David Meyer, Chair 


Although quantum mechanics is linear, there are nevertheless quantum sys¬ 
tems with multiple interacting particles in which the effective evolution of a single 
particle is governed by a nonlinear equation. This includes Bose-Einstein con¬ 
densates, which are governed by the Gross-Pitaevskii equation, which is a cubic 
nonlinear Schrodinger equation with a term proportional to Evolution by 

this equation solves the unstructured search problem in constant time, but at the 
novel expense of increasing the time-measurement precision. Jointly optimizing 
these resources results in an overall scaling of N 1 ^ 4 , which is a significant, but not 
unreasonable, improvement over the N 1 ^ 2 scaling of Grover’s algorithm. Since the 
Gross-Pitaevskii equation effectively approximates the multi-particle Schrodinger 
equation, for which Grover’s algorithm is optimal, our result leads to a quantum 
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information-theoretic bound on the number of particles needed for this approx¬ 
imation to hold, asymptotically. The Gross-Pitaevskii equation is not the only 
nonlinearity of the form f that arises in effective equations for the evo¬ 

lution of real quantum physical systems, however: The cubic-quintic nonlinear 
Schrodinger equation describes light propagation in nonlinear Kerr media with 
defocusing corrections, and the logarithmic nonlinear Schrodinger equation de¬ 
scribes Bose liquids under certain conditions. Analysis of computation with such 
systems yields some surprising results; for example, when time-measurement pre¬ 
cision is included in the resource accounting, searching a “database” when there is 
a single correct answer may be easier than searching when there are multiple cor¬ 
rect answers. These results lead to quantum information-theoretic bounds on the 
physical resources required for these effective nonlinear theories to hold, asymp¬ 
totically. Furthermore, strongly regular graphs, which have no global symmetry, 
are sufficiently complete for quantum search on them to asymptotically behave 
like unstructured search. Certain sufficiently complete graphs retain the improved 
runtime and resource scalings for some nonlinearities, so our scheme for nonlin¬ 
ear, analog quantum computation retains its benefits even when some structure is 
introduced. 
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Chapter 1 


Introduction 

1.1 Linear Quantum Search 

Imagine having a shuffled deck of playing cards where we are searching for 
the Ace of Spades. Since there is no ordering or structure to the cards, one must 
check each card, one by one, until the Ace of Spades is found. It might be the first 
card, or it might be the last card; on average, one must search half of them. If 
there are N cards, then one checks 0(N/2) = 0(N ) of them on average. This is 
the best that a classical computer can do. 

A quantum computer, on the other hand, can solve this problem in 0(y/N) 
steps using Grover’s algorithm [2]. Rather than explaining it in the digital, or 
discrete-time, paradigm in which it was originally proposed, we will focus on its 
equivalent analog, or continuous-time, analogue. This was first given by Farhi and 
Gutmann |3], but we use Childs and Goldstone’s notation and interpretation [I], 

The system evolves in a A r -dimensional Hilbert space with computational 
basis {10),..., \N — 1)}. The initial state |^(0)) is an equal superposition |s) of all 
these basis states: 



N—l 


The goal is to find a particular “marked” basis state, which we label |tu). We do 
this by evolving by Schrodinger’s equation 



1 


2 



Figure 1.1: The complete graph with N — 6 vertices and a single marked vertex 
(colored red). The non-marked vertices are colored white, and the state compo¬ 
nents at them evolve identically by symmetry. 


with Hamiltonian 

Ho = — 7 -/V|s)(s| — |ta)(w|, 

where the first term effects a quantum random walk on the complete graph, so 
7 is a parameter that’s inversely proportional to mass, and the second term is a 
potential well at the marked vertex, causing amplitude to build up there. Since 
the probability amplitudes of finding the randomly walking quantum particle at 


the non-marked vertices evolve identically by symmetry, as shown in figure 1.1 


the system evolves in a two-dimensional subspace spanned by (|ia), |r)}, where 

1 


|r) = 






i^w 


is the equal superposition of the non-marked vertices. 

One might (correctly) reason that the success of the algorithm in finding 
the marked vertex with probability 1 depends on the value of 7 . This can be seen 


in figure 1.2, which shows the difference in eigenvalues of H 0 and the overlaps of 
its eigenvectors with |s) and |iu). When 7 takes a critical value of 7 C = 1/N, the 
Hamiltonian becomes 

H 0 = -|s)(s| - \w)(w\, 


and its eigenstates are 


l^o,i) = —?= 


y/N 


y/2 V VN + 1 


(I S)±\w)) 
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Figure 1.2: Eigenvalue gap and eigenstate overlaps of lit, with N = 1024. 


with corresponding eigenvalues 


E 0 ,i = -1 T 


Vn' 


So the energy gap is A E = 2/y/N. Then the evolution of the system can be 
directly calculated. We begin in the state 


, , . . . 1 y/N + 1 

W°)> = k> = 7i J^ 7r - 


(IV’o) + IV'l))- 


This evolves to 


1 y/N + 1 / 


I m) = e~ iHt \s) = [e~ iE ° t \'il>o) + e~ iElt \i/}i)) 


y/2 V VN 


Plugging in for and |^i), 
1 


I m) = 2 i e ~ lEot (I s ) + k)) + ^ (I s ) - H)} 


= \ [( e~ iEot + e~ iElt 


) I s ) + (' 


j—iEot —iE\t 


) 


w 


= i Le -i(E 0 +E 1 )t/2 ^iAEt/2 + e -iAEt/2 j + ^iAEt/2 _ & -iAEt/2 ) |^j 


= e — h®o+-Ei)^/2 


AE . . . 

cos | —— £ ] |s) + l Sill 


AE 


t I w) 


The amplitude of measuring the randomly walking quantum particle in the vertex 
corresponding to |iu) is 


(«#(*)> = e 


_ p —i{E 0 +Ei)t/2 


COS 


(AE 


_Vn V 2 


t + i sin 


(AE 

\~2 


-t 
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Figure 1.3: Success probability as a function of time for linear search with 7 = 
1/iV. The solid line is N = 100 and the dashed line is IV = 1000. 


So the success probability is 


\(w\^(t))\ 2 = — cos 


/A E 


t + sin 


A E 


which equals 1 when t = n/AE. So the Schrodinger evolution rotates the state 
from |s) to \w) in time n\/N / 2 , as shown in figure 


1.3 


We can also visualize this 


on the Bloch sphere, as shown in figure 1.4, with | w) at the north pole and |r) at 
the south pole; the state starts at |s) near the south pole, moves directly to the 
north pole, loops around the other side, and repeats the motion. 

So the critical 7 is the value of 7 that causes the eigenstates of H 0 to be 
proportional to |s) ± |iu), which causes the system to evolve to the marked basis 
state | w) in 0(\/N) time, thus yielding a successful search. But how do we find 
7 C = 1/N in the first place? Here, we show two methods of finding it. 

The first method of finding the critical 7 is by explicitly finding the eigen¬ 
vectors of Ho and choosing 7 such that they have the desired form proportional to 
|s) ± |u;). Recall the system evolves in the two-dimensional subspace spanned by 
the marked vertex |u>) and the equal superposition of non-marked vertices |r). In 
this basis, the Hamiltonian is 


H 0 = 


-(7 + 1) 

—7-v/iV — 1 


-7v^v^i\ 

-7(^-1)/ 
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Figure 1.4: The evolution of | ip{t)) for linear search with N = 1024 and 7 = 1/iV, 
depicted on a Bloch sphere with |w) at the north pole and |r) at the south pole. 


Let’s find the eigenvalues of H$. The characteristic polynomial is 

det(H 0 - AI) = A 2 + (yiV + 1)A + 7 (N - 1) 

Setting this equal to zero and using the quadratic formula, we get 

—('jN + 1) ± v/( 7 iV - l ) 2 + 4 7 


A = 


which has a gap of 


AA = y/frN - l ) 2 + 47 . 
Now, let’s find the eigenvectors of H 0 : 


IV’i) = 


u 


Writing our coefficient matrix in a more general form, our eigenvalue equation is 

A AA =x f u 

c bj \vj \v, 

which results in (from the second equation) 

A ± — b — 7 N + 27 + 1 =f a/( 7 N — l ) 2 + 47 


u = 


-v = 


27VA - 1 
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or 


' \±—b 


\^±) = 


c 

V 


OC 


'\±-b' 

c 

1 




|iu) + | r) 


A ±-b 


w) + 


N 


N — 1 


s) ~ 


y/N 


\w) 


N 


N — 1 


-'yN + 1 =F y /(7 N - l ) 2 + 47 


oc 


27 s/N 

-'yN + 1 =F \J('yN - l) 2 + 47 

27 \/n 


w) + |s) 


w) + |s). 


We want the term in front of |ta) to equal =)=1. That way, the eigenstates are 
proportional to |s) =F |xu)- This happens when 7 = 7 c = 1/N: 


-7 N + 1 =F \J('yN - l) 2 + 47 


27V A /V 


-T 1 ^ 7 c- N - 


lc 


Another way to determine the critical 7 and runtime is using degenerate 
perturbation theory [3]. We again start with H 0 in the (|w), |r)} basis: 


H 0 = 


-(7 + 1) 
— 'yy/N — 1 


—7 y/N — l\ 


Assuming N is large so that — 1 « N, we separate the Hamiltonian into leading 
order and higher order terms: 


H 0 = 







H, 


(i) 


H, 


( 2 ) 


In lowest order, the eigenstates of the Hamiltonian are |iu) and |r) with correspond¬ 
ing eigenvalues —1 and — 7 N. If the eigenvalues are nondegenerate, then since the 
initial superposition state |s) is approximately |r) for large N, the system will stay 


























7 


near its initial state, never having large projection on \w). For the eigenstates to 
be different, namely a superposition of |r) and |w), we need the eigenvalues to be 
degenerate. That is, when 7 = y c = 1/N, the first-order perturbation causes 
the eigenstates to have the form 

iV’i) = ot w \w) + a r \r ), 


and the coefficients a WjT and eigenvectors E± can be found by solving the eigenvalue 
problem 



where H wr = + H^\r), etc. These terms are easy to calculate. We get 



Solving this eigenvalue problem, we get eigenvectors 


1 

71 



with eigenvalue E + 




with eigenvalue E_ 


Then the approximate eigenstates of E[ 0 are 


1 

Vn 


IV’i) 


1 

71 


(7) =F \r)) 


with eigenvalues 



Note that the energy gap is A E = Since |r) « |s), we approximately have 
the eigenstates from before, so the system evolves from |s) to |tc) in time t* = 
7 t/A E = n\/N /2. 



1.2 Nonlinear Quantum Search 


It is proved that Grover’s algorithm is optimal [5], meaning 0(y/N) is the 
fastest runtime in which quantum mechanics can solve the unstructured search 
problem. To search faster, one must go beyond standard quantum theory, such 
as nonlinear extensions. Abrams and Lloyd [6] gave two examples of nonlinear 
algorithms with fundamental nonlinearities that resulted in unreasonable compu¬ 
tational advantages, solving NP-complete and problems in polynomial time. 
Both of their algorithms can be implemented by a nonlinear Schrodinger-type evo¬ 
lution in which the time derivatives of the state components depend upon their 
hyperbolic tangents BE]. The derivative of tanli x at x = 0 is 1, so this is a 
strongly nonlinear system in which 0 is an unstable fixed point. The strength of 
the nonlinearity provides a large computational advantage, but it also makes the 
system highly susceptible to noise ®nm- 

An obvious question is whether a more modest, physically motivated nonlin¬ 
earity can still produce a computational advantage. While extensive experimental 
work has shown that, at least in the familiar regimes of atomic and optical physics, 
the effect of any fundamental nonlinear generalization of quantum mechanics must 
be tiny PEHE], there are nevertheless quantum mechanical systems with multi¬ 
ple interacting particles in which the effective evolution of a single particle is gov¬ 
erned by a nonlinear equation. These include Bose-Einstein condensates (BECs) 
mmm, whose evolution is described by the celebrated Gross-Pitaevskii equa¬ 
tion [13 [16]: 


.. <9 . , . 

l m ^'( r,t ) = 


h 2 

^ v2+v; 


. 4:Tih 2 a , , , sl2 

extir) -I- N 0 \^(r,t)\ 


m 


V>(r,t). (1.1) 


This nonlinear Schrodinger equation has a cubic nonlinearity, which has zero 
derivative at zero, making it softer than those considered by Abrams and Lloyd. 
In this thesis, we explore the consequences of solving the quantum search problem 
with such a cubic nonlinearity, and we later generalize it to arbitrary nonlinearities 
of the form /(l^l 2 )^, where / is a real-valued function. 

Of course, the Gross-Pitaevskii equation is only an effective approximation 
of the linear, multi-particle dynamics, for which Grover’s algorithm is optimal. So 
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any speedup must be at the expense of increasing the “space” resource such that 
the product of the space requirements and the square of the time requirements is 
lower bounded by N [5], This will yield a lower bound on the number of condensate 
atoms needed for the Gross-Pitaevskii equation to be valid. 

To elucidate the source of the cubic nonlinearity in the Gross-Pitaevskii 
equation, let’s explicitly derive it Ed- The many-body Hamiltonian describing 
multiple interacting particles trapped in an external potential 14 xt (r) with two- 
body interaction potential V(r — r') is 


H = 



h 2 o 

-2^ V +K 


ext ^ 


T(r) + - / drdr / T^(r)T^(r)H(r — r , )T(r / )T(r), 


where we’ve quantized the classical fields by promoting them to creation and an¬ 
nihilation operators, \[b(r) and T(r), respectively (he., second quantization). In 
the Heisenberg interpretation, the state vectors remain fixed while the operators 
evolve according to 




Plugging in H, this becomes 


d 


ih di^^ r,t ^ = 


K 


—V 2 + U ext (r) + / dr , T t (r,t)H(r - r')T(r / ,t) 


T(r,t). 

We express the operator using mean held theory as an order parameter plus a 
perturbation: 

T(r,t) = (T(r, t)) + T'(r ,t). 


The order parameter, or number density, can be normalized and interpreted as 
the wave function of the condensate, so we write it as -0(r ,t) = {^(r,t))/y/No, 
where N 0 is the number of condensate atoms. Assuming that the perturbation 
is negligible, so the temperature of the condensate is near 0, we get T(r ,t) —)• 
'ijj{r,t)y/No. When the Bose gas is dilute, meaning the s-wave scattering length 
a is much less than the interparticle spacing, then the effective interaction is (see 
section 5.2.1 of [IS] ) 


t/y o 4vrh 2 a 
V (r — r ) = -o(r — r ). 


m 


Using this, the evolution of the condensate wave function becomes 


.. d 

lh 0l U( r, G = 


h 2 9 

~^ v+v ‘ 


ext 


r + 


4vrh 2 c 


m 


■No |^(r,t)| J 


^( r ,t), 
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which is the Gross-Pitaevskii equation in Eq. |1.1 

Physically, we can exert some control over the strength of the cubic non¬ 
linearity in the Gross-Pitaevskii equation by varying the scattering length a via 
Feshbach resonance [15] , In this process, two condensate atoms interact via the 
hyperhne interaction (he., an interaction between the electronic and nuclear spins 
of the atoms), forming a quasi-bound state. The energy of this unstable interme¬ 
diate state is higher than when the atoms are separate by the binding energy Ey n 
and it can be further offset with a “detuning” e. When e = 0, it is said that the 
collision is “on resonance.” In an optically trapped BEG, the detuning corresponds 
to an external magnetic field B [20], and the effective scattering length near the 
resonance Bq is 



where a is the scattering length away from the resonance and A B is related to 
the width of the resonance [19]. Thus there is no theoretical limit as to how much 
the scattering length can be varied using Feshbach resonance. Experimentally, it 
depends on the precision in which the external magnetic field can be controlled, 
and the first group to experimentally observe Feshbach resonance in BECs was 
able to vary the scattering length by a factor of 10 m■ 

The first BEC to be experimentally produced was made by a team led by 
Eric Cornell and Carl Wieman of JILA by trapping and cooling 2 x 10 4 rubidium- 
87 atoms in a magnetically-confined trap |2lj. About four months later, this was 
improved by Wolfgang Ketterle’s team at MIT, who trapped 5 x 10 5 sodium atoms 
with the addition of an optical plug [22]. Both rubidium-87 and sodium atoms have 
positive scattering lengths, meaning the bosons are repulsive. While BECs with 
negative scattering lengths cannot exist as a homogeneous gas since condensation 
is preempted by a first-order phase transition [23], they are stable against collapse 
in the non-homogeneous environment of a trap for small numbers of atoms less 
than a critical value given by 


a 


where k is the dimensionless “stability coefficient” depending on the ratio of mag¬ 
netic trap frequencies, and Oho is the harmonic oscillator length [2TI [25]. This was 
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experimentally verified when a condensate of roughly 10 3 lithium-7 atoms was pro¬ 
duced by Randy Hulet’s team at Rice University, just one month after the JILA 
collaboration’s discovery [26J E2!. Condensation of bosons with attractive interac¬ 
tions can also be demonstrated using Feshbach resonance, using the detuning to 
turn the interaction from repulsive to attractive [ 28} . 

So the Gross-Pitaevskii equation is rooted in established physics, mak¬ 
ing its cubic nonlinearity a physically reasonable term to include in computa¬ 
tion. In the next chapter, we quantify the computational advantage that the 
cubic nonlinearity provides for the unstructured search problem compared to stan¬ 
dard quantum computation. This requires introducing a novel physical resource: 
time-measurement precision. Since this advantage cannot persist when the Gross- 
Pitaevskii equation is recognized as an approximation to an underlying multi¬ 
particle Schrodinger equation, for which Grover’s algorithm is optimal, we arrive 
at a quantum information-theoretic lower bound on the number of condensate 
atoms needed for this approximation to hold, asymptotically. 

In Chapter 3, we generalize nonlinear search on the complete graph to 
arbitrary Schrodinger-type nonlinearities of the form /(|^|-’|^), where / is a real¬ 
valued function. This includes the cubic nonlinearity in the Gross-Pitaevskii equa¬ 
tion, as well as other physical systems including the cubic-quintic and logarithmic 
Schrodinger equations. This yields some surprising results; for example, when 
time-measurement precision is included in the resource accounting, searching a 
“database” when there is a single correct answer may be easier than searching 
when there are multiple correct answers. 

As previously explained, search on the complete graph evolves in a two- 
dimensional subspace spanned by the marked vertex |u>) and the superposition of 
non-marked vertices |r), which we colored red and white in figure 0 The next 
level of difficulty is search in a three-dimensional subspace, which strongly regular 
graphs support. In Chapter 4, we use degenerate perturbation theory in a novel 
way to solve the quantum search problem on strongly regular graphs, showing that 
search also achieves the 0(\/ r N) speedup. This is similar to the hypercube, which 
evolves in a larger space than the complete graph, but still searches in 0(yfN) time 
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PP. Search on these “sufficiently complete” graphs is sped up by the nonlinearities 
in the same way that that it is for the complete graph in Chapters 2 and 3, and 
we work through this in Chapter 5. 

Finally, we conclude with a summary and give some future directions in 
Chapter 6. 


Chapter 2 


Nonlinear Quantum Search on the Complete 
Graph 

2.1 Setup 

To review the search problem, the system evolves in a TV-dimensional Hilbert 
space with computational basis {|0),..., \N — 1)}. The initial state |^(0)) is an 
equal superposition |s) of all these basis states: 

i N - 1 

\m) = w = ^ £ k>. 
v 2=0 

The goal is to hnd a particular “marked” basis state, which we label |tc). 

In the nonlinear regime, we include an additional nonlinear “self-potential” 
V(t) so that the system evolves according to the Gross-Pitaevskii equation Eq. m 

[H 0 - g\i/j(r : t)\ 2 ]i/j(r,t), 

where g > 0. This corresponds to a BEG with attractive interactions, and thus 
a negative scattering length [S3J |2E] • Heuristically, as probability accumulates at 
the marked state due to the |u')(w| term in Hq, the self-potential attracts more 
probability, speeding up the search. Thus we expect larger g to result in a faster 
algorithm. 

In the computational basis, the self-potential is 

N -1 

V(t) = 9^2 l(#)l 2 N>(*l- 

»=o 
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Figure 2.1: Success probability as a function of time for N = 1024 and 7 = 1/N 
constant. The solid curve is the linear (g = 0) case, and the dashed curve is the 
nonlinear g — 1 case. 


Even with this nonlinearity, the system remains in the subspace spanned by (|ta), 
\s)} throughout its evolution. We define a vector 


r) 


1 




which is orthonormal to \w). Then the state of the system \ift(t)) can be written 
as 


IV’W) = a(t)\w) +P(t)\r). 


Writing the Gross-Pitaevskii equation in this {|w), |r)} basis, we get 
d I a 


a 


= -1 (H 0 - V ) 
dt \pJ \P. 


= 1 


. (1 + 1 + g\a\ 


7 \/N — 1 


a 


7^V^1 7 (iV-l) + j^l/JlV W’ 


( 2 . 1 ) 


A 


where we’ve defined A = —(H 0 — V). 


2.2 Critical Gamma 

Before proceeding with further analytical calculations, we build some intu¬ 
ition by examining two plots. For constant 7 and g, the success probability as a 
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Figure 2.2: The evolution of \'i/(t)) for nonlinear search with N = 1024, 7 = 1/N 
constant, and g — 1, depicted on a Bloch sphere with \w) at the north pole and 
|r) at the south pole. The system fails to reach the north pole, and so it fails to 
reach a success probability of 1. 



Figure 2.3: Eigenvalue gap and eigenstate overlaps of A at t = 20 for nonlinear 
search with N = 1024, g = 1, and 7 = 1/N constant. 
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function of time, \a(t)\ 2 , is plotted in figure 2.1 along with the linear result. The 
nonlinear algorithm underperforms the linear one in this case. As shown on the 


Bloch sphere in figure 2.2, with | w) at the north pole and |r) at the south pole, 
the system starts near the south pole and begins moving towards the north pole. 
But then it veers to the side, looping near the bottom of the sphere, and returning 
to its initial position. So the system never has high success probability. This is 
true in general for constant 7 and g, and it can be understood by examining the 
time-dependence of the critical value of 7, which is the value of 7 that ensures 
that the eigenstates of A are in the form ±|ta) + |s). Initially, y c = 1/N. Then, as 
shown in figure [273 , it shifts to a larger value. If 7 is constant, it will not follow 
this shift, we will no longer have the desired eigenstates, and the algorithm will 
perform poorly. 

To determine how y c varies with time, we find the eigenvectors of A and 
choose 7 so that they have the desired form ±|w) + |s). To eliminate fractions in 
the subsequent algebra, we rescale the nonlinearity coefficient g by defining 

r- 9 

Solving the characteristic equation gives the eigenvalues of A: 

A± = ^ ( 7 N + 1 + Go) ± ^AA, 

where the gap between them is 

AA = v/( 7 A - l ) 2 + 47 + G 2 S 2 + 2 GS [1 - 7 (IV - 2)], 
and we’ve defined 

o = (N — l)\a\ 2 + \/3\ 2 and 5 = (N — l)|a | 2 — \/3\ 2 . 


The corresponding eigenvectors of A are 

7 IV + 1 + <5G ± AA 


IVa) = 


N 

N — 1 


27 \/N 


w) + | s) 


The critical value of 7 ensures that these eigenvectors have the form ±|w) + |s). 
That is, 


— 7 N + 1 + 5G ± AA 


2 r y\f~N 


= 1 . 


7c 
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Figure 2.4: Success probability as a function of time for nonlinear search with 
G — 1 and 7 = y c as defined in Eq. IT2 The solid line is IV = 100 and the dashed 
line is N = 1000. 


Solving this yields: 


7c = 


1 + G5 
N ' 


( 2 . 2 ) 


Note that in the linear limit (G — 0), this reduces to y c = 1/N, as expected and 


calculated in Chapter 1. Since <5 varies with time, Eq. T2 implies q c also varies 
with time, in agreement with our previous discussion about figures |2.1 and 2.3 


Furthermore, it can be precomputed without needing to know the location of the 
marked vertex, so there is no issue of having to measure the system during the 
computation. 


2.3 Runtime 

For the remainder of the chapter, we choose time-varying 7 = q c according 
to Eq. B Before analytically determining the consequences of this, let’s again 
consider a plot. Figure [2N shows the success probability as a function of time. 
There are several observations. First, the success probability reaches 1, which 
occurs because we constructed the eigenstates to make this happen. Second, as N 
increases, the runtime remains constant. Third, the success probability is periodic. 
Finally, the peak in success probability becomes increasingly narrow for large N. 
Let’s now analytically prove the second, third, and fourth observations. 
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To begin, we explicitly write out Eq. |2.1| to get two coupled, first-order 
ordinary differential equations for a(t) and (3(t): 


^ = i {[7c + 1 + G(JV - l)|a| 2 ] a + 7c Vivm/s} (2.3) 

^ = * {y c VW=la + [ 7c (JV - 1) + G\t3\ 2 } p) . (2.4) 


We decouple these equations by defining three real variables x(t), y(t), and z(t ) 
such that 


x = \a\ 
y + iz = a/3*. 


(2.5) 

( 2 . 6 ) 


Note that x(t) defined by Eq. 2.5 is the success probability. Differentiating it and 


utilizing Eq. 2.3, we find that 

dx d\a\ 2 

dt dt 

Solving this for z, we get 


da* da . r ——- 

= a— r— H ——a = 27 c \/ N — 1 z. 


dt 


dt 

1 dx 


z = 


27 c\/lV - 1 dt' 


(2.7) 


Noting that d^Jdt = Gdx/dt, we differentiate Eq. 2.7 to get 

2 1 d 2 x 


dz 

dt 


2 v / F ^1 7, 2 


I + 


7 c dt 2 


( 2 . 8 ) 


Now we want to find another expression for dz/dt , which we can then set equal to 
Eq. [275} We do this by differentiating Eq. |2.6[ utilizing Eq. |2.3| and Eq. |2.4[ and 
equating the real and imaginary parts, which yields 


d V o 

dt = “ 27 - 

= ‘ilcV + ~tcVN - 1(1 - 2j;j 


(2.9) 

( 2 . 10 ) 


Substituting Eq. 2.7 for z into Eq. 2.9, we get 


dy —1 dx 
dt y/N — 1 dt 1 
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which integrates to 


V = 


0 -*): 


where the constant of integration was found using x(0) = 1/N and y(0) = 
yj N — 1/N. Now we can plug this into Eq. 

1 


2.10 


to get 


dz 

dt ~ T 


(1 — x) + jcVN — 1(1 — 2x ) 


7c 


y/N^l 


(1 + N — 2Nx) 


Equating this to Eq. 2.8 and simplifying yields 




Plugging in for y c as defined in Eq. 2.2, this becomes 

W = 1 -G + NGx(^) +^ 2 ( 1 - g + ng *) 2 ( 1 + n ~ 2n *) P- 11 ) 


Now let f(x) = ( dx/dt ) 2 so that df/dx = 2d 2 x/dt 2 . Then Eq 

1 d[_ NG x , 2 / 1 ^ 

2 dx 


2.11 


becomes 


NG 2 

Y -C Tmrx f + W 2 (l-G + NG X r(l + N-2N X ). 


Solving this first-order ODE and using the initial condition f(x = 1/N) = 0, we 
get 


fix) = 


4(iVx - 1)(1 - x) [1 + G(Nx - 1)]' 
N~ 2 


Taking the square root and noting that dx/dt = ±y/f(x), 


dx 

Tt =± ’ 


U{Nx - 1)(1 - x) [1 + G(Nx - 1 )f 
N 2 


( 2 . 12 ) 


To solve this uncoupled equation, we use separation of variables and integrate from 
t — 0 to t and x = 1/N to x, which yields 

y/N yjl — x 


t = 


N 

1 + G{N-1) 


tan 


-l 


y/l + G{N -l)y/Nx-l 


TT 

2 


(2.13) 


Solving for x, the success probability as a function of time is 

N + [1 + G(N - 1)] tan 2 


x{t) = 


7T l+G(jV~l) . 

2 V N L 


N + N [1 + G(N - 1)] tan 2 


l+G(Af—1) 


N 


(2.14) 
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From this, the success probability reaches 1 when the tangent term is zero, which 
first occurs at time 




7r 


s/N 


y/l + G(N-l) 2 

This runtime is exactly constant for G = 1. Also, when G = ©(1), the runtime 
for large N is tt/2\/G, and thus asymptotically constant (and arbitrarily small!). 
From Eq. 2.14[ we also see that the success probability is periodic with a period 


of 2 f*. 


Now let’s prove that the peak in success probability is narrow by finding 


its width, thus proving all our observations about figure 2.4 Using Eq. 2.13 the 
difference in time at which the success probability reaches a height of 1 — e is 


At 


/ N tan-' 

VNyfe 

/ 1 + G(N — 1) 

y/1 + G(N — l)y/N(l - e) — 1 


The tan 1 makes it difficult to determine the scaling with N, so we Taylor expand 
it: 

Ai = 1 + g{N -1) VftT + 0(£3/2) ' 

When G = N K , the first term scales as ©(IV 1 / 2 ) when k < — 1 and ©(iV -1 / 2 ^) 
when k > — 1, for large N. To determine whether keeping this first term alone is 
sufficient, we use Taylor’s remainder theorem to bound the error 


Ri(e) < 


N 2 (l + 3G(N(1 — e) - 1)) 


(N( 1 - e) - 1 ) 3/2 (1 + G (N (1 - e) - l)) s 


U/2 


which has the same scaling for large N as the first term in the Taylor series for 
At. Thus it suffices to keep only the first term. 

For constant G , the width in success probability is ©(l/x/iV), which agrees 
with our observation from figure |2~4| that the peak in success probability is increas¬ 
ingly narrow as N increases. Thus we must measure the system with increasing 
time precision. This behavior is opposite the linear case. That is, when G = 0 the 
width is &(y/N), so the time at which we measure the result can be increasingly 
imprecise as N increases. 
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2.4 Time-Measurement Precision 


This time-measurement precision requirement of the nonlinear algorithm 
requires additional resources. In particular, time and frequency standards are cur¬ 
rently defined by atomic clocks, such as NIST-F1 in the United States [29] . An 
atomic clock with n c i oc k ions used as atomic oscillators has a time-measurement 
precision of 1 / y/n c i oc k when the ions are acted upon independently. This can be 
improved using quantum entanglement, reducing the time-measurement precision 
to l/n c i oc k [3Di GSj. Even with this improvement, our constant-time nonlinear 
search algorithm would require 0(\/TV) ions in an atomic clock to have sufficiently 
high time-measurement precision to measure the peak in success probability. So, 
although our nonlinear algorithm runs in constant time, the total resource re¬ 
quirement is still 0(VN), the same as the linear algorithm. This raises the pos¬ 
sibility that nonlinear quantum mechanics may not provide efficient solutions to 
NP-complete and problems when all the resource requirements are taken into 
consideration [6j. 

In our case, however, we can settle for a smaller improvement in runtime 
and reduce the time-measurement precision and total resource requirement. If we 
let G decrease as N K for k < 0, then the runtime is £* = @(TV _K / 2 ), and the 
time-measurement precision is At = @(TV~ 1 / 2 ~'' i ), where we’ve assumed for both 
that k > —1, since for k < —1, At = ©(iV 1 / 2 ), independently of G. This time- 
measurement precision requires 0(N 1 ^ 2+K ) ions in an atomic clock that utilizes 
entanglement. We assume, as in the setup for Grover’s algorithm, that log TV 
qubits can be used to encode the TV-dimensional Hilbert space; these should also 
be included in the required “space” resources. Multiplying the time and “space” 
requirements together, which preserves the time-space tradeoff inherent in naive 
parallelization, the resulting total resource requirement takes a minimum value 
of O (TV 1 / 4 log TV) when k = —1/2 (so that the runtime is TV 1 / 4 and the time- 
measurement precision is constant). The success probability as a function of time 


at this jointly optimized value of G is plotted in figure 2.5; note that the peak 
width is independent of TV. 

This significant—but not unreasonable—improvement over the 
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Time 


Figure 2.5: Success probability as a fu nctio n of time for nonlinear search with 
G = iV -1 / 2 and 7 = y c as defined in Eq. 2.2 The solid line is N — 100 and the 


dashed line is IV = 1000. The peaks have same width, independent of N. 


O(VNlogN) time-space resource requirements of the linear quantum search algo¬ 
rithm is consistent with our expectation that a modest nonlinearity should result 
in a modest speedup. 


2.5 Repulsive Interactions 

Our nonlinear search algorithm was based on the intuition that attractive 
interactions speed up the accumulation of success probability. By the same in¬ 
tuition, repulsive interactions, where G < 0, should yield a worse runtime. Our 
derivation of Eq. |2.2| for the critical value of 7 is unchanged if we flip the sign of 


G, so Eq. 2.11 and Eq. |2.12 are still valid for repulsive interactions. These equa¬ 
tions yield critical points x * = 1 /N, 1, and (G — 1 )/NG, corresponding to minima, 
maxima, and stationary points, respectively. 

When G > —1/(N — 1), the success probability is unhindered by the sta¬ 
tionary point and reaches a maximum value of 1 , as shown in the dashed curve of 


figure 2.6. When G < —1/(N — 1), however, reaching this maximum is precluded 


by the presence of a stationary point, as shown in the dashed and dot-dashed 


curves of figure 2.6 


We can explicitly prove that repulsive interactions (G < 0) will underper- 
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Figure 2.6: Success probability as a function of time for N = 1024 and 7 = y c 

The solid line is the linear (g = 0) case, the dashed line is 
0.5 case, the dotted line is the nonlinear g = —1 case, and the 


as defined in Eq. 2.2 
the nonlinear g = 


dot-dashed line is the nonlinear g = —1.5 case. 


form the linear (G = 0) algorithm. From Eq. 2.12 


d 'T' 0 

— = ±-^(Nx- !)(! -x) [1 + G(Nx 


!)]• 


So when G < 0, the magnitude of dx/dt at a particular value of x is less than when 
(7 = 0. Then success probability will increase more slowly for repulsive interactions 
(except initially, where they increase at the same rate). Thus it will underperform 
the linear algorithm. 


2.6 Validity of the Gross-Pitaevskii Equation 

Of course, the cubic nonlinearity we’ve exploited is not fundamental, but 
rather occurs in an effective description of an interacting multi-particle quantum 
system (e.g., a BEC). So we must include the number of particles Nq in our resource 
accounting. Each particle interacts with the potential at the marked site, so in 
the framework of Zalka’s optimality proof for Grover’s algorithm [5j (generalized 
to continuous time [32j), there are N 0 oracles, each responding to a logiV bit 
query. Zalka showed that the product of the space requirements and the square 
of the time requirements is lower bounded by N, i.e ., (No logiV)(iV 1 / 4 ) 2 = Q(N). 
Solving for the number of particles, Nq = 0(iV ly/2 / log N). This is a quantum 
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information-theoretic lower bound on the number of particles necessary for the 
Gross-Pitaevskii equation to be the correct asymptotic description of the multi- 
particle (linear) quantum dynamics. 

Notice that once we account for the scaling of N 0 in the space require¬ 
ments, the product of the time and space requirements is 0(iV 3 / 4 ) ; worse than the 
0(N 1 ^ 2 log TV) of Grover’s algorithm. In fact, if we calculate for the general case 
G = N K , where k need not be chosen to optimize the product of the time and space 
(ignoring N 0 ) resources, Zalka’s bound implies N 0 = G(max{l, N 1+K / logiV}), so 
the total time-space requirements are 0(N 1+K / 2 ) for k > —1, and 0(N 1 ^ 2 log N) 
when k = —1. This is optimized for k = —1, he., by Grover’s algorithm. On the 
other hand, Zalka’s bound is strongest when k — 0, in which case it implies that 
Nq = Vt(N/ log TV). That is, the existence of the constant time nonlinear algo¬ 
rithm we found in section 4 implies this stronger lower bound on N 0 , despite the 
0(N x / 2 ) number of clock ions required. To our knowledge, this is the first lower 
bound derived on the scaling of N 0 required for the Gross-Pitaevskii equation be 
a good asymptotic approximation. 

This bound also is significantly stronger than the bound implied by the 
physically plausible requirement that the volume of the multi-particle condensate, 
and thus N 0l be of at least the order of the volume of space in which the N possible 
discrete locations are defined. Were we working in any fixed, finite dimension, e.g., 
on a cubic lattice, the volume would be proportional to N, implying N 0 = f 1(N). 
But we are not; the complete graph with equal pairwise transition rates is realized 
by the vertices and edges of an equilateral (N — 1)-dimensional simplex. With 
edges of length 1, this has volume \J iV/2 Ar_1 / (N — 1)!, which is much smaller than 
N, and also much smaller than our bound of N/ log N. 

2.7 Critical Gamma is a Continuous Rescaling of Time 

We previously derived the critical value of 7 so that the eigenstates of the 
Hamiltonian are proportional to ±|w) + |s). Now we examine what the critical 
value of 7 does from another perspective. Recall the “Hamiltonian” we’ve been 
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using is 


H = — 7 _/V|s)(s| - \w)(w\ \ipi\ 2 \i){i\, 


where ipi = Explicitly writing the nonlinear term as marked and unmarked 

terms, we get 


| «i 2 

H = — y.V .s)(.s] - \w)(w\ - g\a\ 2 \w)(w\ - gjp~^ l x X 

X^W 


X\ 


— 7 iV|s)(s| — \w)(w\ — G(N — l)|a| 2 |ta)(w| — G\/3\ 2 ^ \x)(x\ 


X^W 


- 7 -ZV|s)(s| — [l + G(N — l)|a| 2 ] \w){w\ — G\f3\ 2 ^ |x)(x|. 

X^W 


Recall 7 = 7 C is chosen according to Eq. 2.2 


7 c-^ = 1 + G(N — l)|a | 2 — G\/3\ 2 , 


which we arrange to get 


1 + G(N - l)\a\ 2 = >y c N + G\p\' 


Then the Hamiltonian becomes 


H = -7ciV|5)<5| - [ lc N + G\p\ 2 ] \w){w\ - G |/3 | 2 |x)(x| 

X^W 

= - 7 f c N (|s)(s| + \w)(w\) - G\f3\ 2 \. 


The last term continuously redefines the “zero” of energy, so we can drop it. That 
is, it only changes the overall phase of the system, which has no measurable effect. 
Then the Hamiltonian is 


H = — 7 IV (|s)(s| + |w)(w|). 

Importantly, H^q = — |s)(s| — |ia)(m| is the Hamiltonian from Farhi and Gutmann’s 
“analog analogue” of Grover’s algorithm J3], and it is optimal. Our nonlinear 
algorithm has a factor of 7 N, so it effectively follows their optimal algorithm, but 
with a continuously rescaled time. That is, the system evolves according to 

. dtp 
'yNdt 


Hfg 
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Let’s call the rescaled time r(t) so that dr = jNdt. Then 


r = / 7 Ndt, 


and the equation of motion becomes 

M rr , 

i— = HpQlp. 

dr 

This has success probability given by (11) of |3j: 


x(t) = sin 2 


r 


\ 1 

+ TT cos 


y/N J N WnJ' 


Plugging in for r, 


x(t) = sin 2 


2 ( 1 2 ( jl N dt\ 


- + — cos - 


. 


V VN J N V y/N ) 

Since 7 C N = 1 + G5 = 1 — G + GNx , we get 

’ (1 — G)t + GN f x(t)dt 


. 2 f 0- — G)t + GN Jx(t)dt^ 1 2 f( 

x(t) = sm (- 7n - J + ^ cos (' 


\/N 


This integral transcendental equation gives x(t). While the difficulty of solving this 
equation makes it less useful in practice, it does reveal our nonlinear algorithm’s 
relationship with the linear, optimal algorithm. In particular, a different control 
policy for 7 will cause the system to evolve along a different, slower path. While 
not a proof, this is an argument for the optimality of our algorithm. 


2.8 Multiple Marked States 


Our analysis naturally extends to the case of k marked states. Let M be 
the set of marked basis states. As before, the system evolves in a two-dimensional 
subspace: 


lV>(t)) = a(*)4= + /J— 

V k ttM - 


i x )- 

x(£M 


The system evolves according to 
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where 

A = 7iV|s)(s| + U + 9^j~] Y1 I^X^I 

' ' x eM x£M 

includes both the linear Hamiltonian and the nonlinear “self-potential”. The eigen¬ 
states of A have the form ±|w) + |s) when 7 is 

1 + G5 


where G = g/(k(N — k)) and 5 = (N — k)\a\ 2 — k\/3\ 2 . At 7 = j c , we can decouple 
these equations in the same manner as the k — 1 case and integrate from t — 0 to 
t and x = k/N to x to get 


t = 


N 


k(l + G(N — k)) 


tan 


-1 


VNy/l 


X 


\Jl + G(N — k)y/Nx — k 


which can be solved for a success probability of 


N T k [1 T G(N — k)] tan 2 

7T 

2 A 

/ k(l+G(N-h)) t 

N + N [1 + G(N - k)] tan 2 

7T 

2 

/k(l+G(N-k)) t 


Then the runtime is 

+ 1 ny/N 

U ~ y/k(l + G(N 


and the success probability is still periodic with period 2f*. At this runtime, the 
peak in success probability has a width of 


1 N -1 

VNyfe 

]k(l + G(N-k)) an 

y/l + G(N - k)y/N( 1 -e)-k 


but Taylor’s theorem can be used to show that it suffices to keep the first term in 
the Taylor series: 

At = 1 + G(N — k) ]/k(N-k) + °( e3/2 )‘ 

As in the case of a single marked state, we can find the scaling of G — N K 
that optimizes the product of “space” and time, where “space” includes both the 
number of ions needed in an atomic clock that utilizes entanglement to achieve 



sufficiently high time-measurement precision, and the log IV qubits needed to en¬ 
code the TV-dimensional Hilbert space. Say the number of marked sites scales as 
k = N x , with 0 < A < 1. When k = —A/2 — 1/2, the product of “space” and 
time takes a minimum value of ST = jV _A / 4+1 / 4 log N (so that the runtime is 
jy— A/4+1/4 an( ] |-q ie time-measurement precision is constant). Note this is a square 
root speedup over the linear (G = 0) algorithm, whose product of “space” and 
time is /v _A / 2+1 / 2 logTV. Thus our nonlinear method, by varying 7 and choosing 
an optimal nonlinear coefficient G, provides a significant, but not unreasonable, 
improvement over the continuous-time analogue of Grover’s algorithm, even with 
multiple marked items. 

Chapter 2, nearly in full, is a reprint of the material as it appears in “Non¬ 
linear Quantum Search Using the Gross-Pitaevskii Equation” in New Journal of 
Physics 15, 063014 (2013). D. A. Meyer and T. G. Wong both contributed signif¬ 
icantly to the work. 



Chapter 3 


Quantum Search with General Nonlinearities 

3.1 Introduction 

So far in this thesis, we’ve only considered Schrodinger evolution with a 
cubic nonlinearity, i.e., evolution by the Gross-Pitaevskii equation [13 fib] : 



where H 0 includes the kinetic energy and trapping potential, m is the mass of the 
condensate atom, and Nq is the number of condensate atoms. In Chapter 2, we 
quantified the computational advantage that this cubic nonlinear Schrodinger equa¬ 
tion has in solving the unstructured search problem. To summarize, we search for 
one of k “marked” basis states among N orthonormal basis states {|0), 11),..., \N — 
1)}. Without the nonlinearity, the optimal solution is the continuous-time ana¬ 
logue of Grover’s algorithm HElElEi, which runs in time N/h). With the 
nonlinearity, we can search in constant time with appropriate choice of parame¬ 
ters, as shown in figure |2.4 This figure also reveals that the success probability 
spikes suddenly, so increasingly precise time measurement is necessary to catch 
the spike. This requires a certain number of atoms in an atomic clock that utilizes 
entanglement ]311 1 .30] . Jointly optimizing the runtime and number of clock ions, 
we achieve a resource requirement of 0((N/k ) 1 ^ 4 )—a square-root speedup over the 
linear quantum algorithm. 

As explained in Chapter 2, Grover’s algorithm is optimal [5], so there must 
be additional resources such that the product of the space requirements and the 


29 
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square of the time requirements is lower bounded by N. In the case of the Gross- 
Pitaevskii equation, the additional resource is the condensate atoms, and the bound 
on the number of them is strongest at Q(N/logN) for the constant-runtime al¬ 
gorithm. Thus we’ve found a quantum information-theoretic lower bound on the 
number of condensate atoms needed for the Gross-Pitaevskii equation to be a good 
asymptotic description of the many-body, linear dynamics. 

These two results—a significant, but not unreasonable, square-root speedup 
in solving the unstructured search problem, and the lower bound on the resources 
necessary for the Gross-Pitaevskii equation to be valid—suggest it is valuable to 
quantify the computational advantage that other effective nonlinearities have in 
solving the unstructured search problem. In particular, we consider nonlinear 
Schrodinger equations of the form 

i w = i H °~ 9 ^ i^i 2 )]^’ l 3 - 1 ) 


where / is some real-valued function. The cubic nonlinear Schrodinger equation is 
the case when f(p) = p. 

A reasonable way to adjust the cubic nonlinearity is to include higher-order 
terms, such as the quintic term that appears when three-body interactions are 
included in the description of a BEG [33j. Another example of including higher- 
order terms is the propagation of light in Kerr media |3T, 35] [36], whose quantum 
origins are worked out in [37] • When a material is subjected to an electric field E, 
its index of refraction n changes: 


, dn 1 d 2 n 


But from symmetry, many materials require that the index of refraction be an even 
function. Then the first-order term is zero, leaving 

. 1 d 2 n 9 

n W= n+ 2d& E + '~ 


The electric field needn’t come from an external source—it can be from the incident 
light itself. For certain incident light beams, this second-order correction is self- 
focusing, and it appears in the equation of motion as a cubic tern[|] The cubic 

1 Since the intensity is proportional to the square of the electric field, the index of refraction 
is frequently written as n(I) = no + 712/. 
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self-focusing term, however, is sometimes insufficient to describe the propagation, 
and a quintic defocusing correction must be included [3S113U1 HO] . This results in 
the cubic-quintic nonlinear Schrodinger equation 

i ~J+~ = -l N ^n ~ g (| ^n \ 2 3 ~ IV'nl 4 ) Ipn, 


which naturally describes a periodic array of N waveguides, where 7 is a parameter, 
ip n is the amplitude of the electromagnetic wave in each waveguide, and A is 
the discrete second derivative [31]. This equation is of the form of (3.1) with 
fip) =P~P 2 - 


The above nonlinearities, and indeed general nonlinearities of the form (3.1), 
do not retain the separability of noninteracting subsystems. That is, in (linear) 
quantum mechanics, if a physical system consists of two noninteracting subsys¬ 
tems, then its state can be written as the product of the states of the subsystems 
(he., as a product state). Nonlinearities, however, generally cause initially uncor¬ 
related subsystems to become correlated. The one exception [32j is the special case 
when f{p) = log(p). Then separability is retained, and the nonlinear Schrodinger 


equation (3.1) contains a loglinear term: 


^ = [-Ho -#log(|^| 2 )] ip - 


Note that the limit of y/x log(a;) as x goes to 0 is 0, so the evolution doesn’t cause 
the wavefunction to diverge^ Not only is the logarithmic^ nonlinear Schrodinger 
equation important for its uniqueness in retaining separability, but it may be suit¬ 
able for describing Bose liquids, which have higher densities than BECs [33]. 

In the following section, we write the generalized nonlinear search problem 
with multiple marked vertices in its two-dimensional subspace. Then we solve it, 
referencing the solution to the cubic nonlinear Schrodinger equation from Chapter 
2 as we go. Finally, we end with two comprehensive examples of searching with 
the cubic-quintic and loglinear nonlinearities that were introduced above and give 
lower bounds for the physical resources needed for them hold. 

2 This concern was also addressed in l l.'ll by examining the generalized Lagrangian density and 
effective potential density in [ 32 ] . 

3 Although the nonlinearity is loglinear, the equation is typically referred to as logarithmic. 

This is different from the cubic and cubic-quintic nonlinearites where the equations are also 
referred to as cubic and cubic-quintic, respectively. 







32 


3.2 Setup 


On the complete graph of N vertices, we have k marked vertices, of which 
we are looking for any one. Let’s call the set of marked vertices M. Then the 


nonlinear Schrodinger equation (3.1) has 


H 0 = —7-/V|s)(s| - Y 


ieM 


from which we subtract a nonlinear “self-potential” 


N—l 


v (t) = f/^/(l(#)| 2 ) l*)(*l 


i =0 

As with the cubic nonlinearity in Chapter 2, g must be positive for the nonlinear 
algorithm to perform better since, heuristically, it causes the self-potential to act 
as an additional potential well, therefore attracting more probability and speeding 
up the search. 

As the system evolves, it remains in the two-dimensional subspace spanned 
by orthonormal vectors 


1 

a Jk 


|i) and 






ieM v " i£M 

so we can write \if(t)) as a linear combination of them: 

1 ^ _ 1 


m)) = a[t) v% 5 , M + m 7W=% 


x£M v ^ x$M 

Then the probability of measuring the system in basis state \i) is 


l <#>| 2 = 


M 2 

k ’ 

i G M 

l/?l 2 

JV-fc’ 

i <£ M 


Let’s define 


fa — f 


|a| 

~k 


and f /3 — f 




N-k 


Then the nonlinear Schrodinger equation (3.1) is written in the two-dimensional 
subspace as 


d (a 


a 


= i 


. f'yk + 1 + gf a 'fVky/N - k 
yys/ky/N - k ry(N - k) + gfpj \/3 


(3.2) 
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3.3 Critical Gamma 


We can also write H{t) = H 0 — V(t) in terms of f a and fp\ 

h = -7iv|s)(s| - Y !*)(*! -gf<*Y !*)(*! -9UYm\ 


iGM 


i&M 


i£M 

N—l 


= -7^V|s)(s| - (1 + gf a - ~ 9f?Y !*)(*!■ 

*SM i=0 

The last term is a multiple of the identity matrix, which simply redefines the zero 
of energy (or contributes an overall, non-observable phase), so we can drop it. 
From our previous work on the cubic nonlinear Schrodinger equation in Chapter 
2 , the critical 7 causes the nonlinear system to follow the same evolution as the 
linear, optimal algorithm, but with rescaled time. That is, we choose 


7 c = 71. [1 + g (/«- M ] = [1 + 9 (/<» - M \ , 


(3.3) 


which is time-dependent, so that 


H = (1 + gfa - gfp) ( ~ 1 lN\s)(s\ - Y l*)(*l ) = (! + gf<* - gfp) #fg- 


i£M 


So the system evolves according to Farhi and Gutmann’s Hamiltonian, but with 


continuously rescaled time. Thus we have the critical 7 (3.3) for general nonlin¬ 


earities of the form ( 3.1[ ). Note that for the cubic nonlinearity, f(p) = p, so if we 
define G = g/[k(N — k)] and 5 = (N — k)\a\- — k\/3 \ 2 , then we get the familiar 


result (1 + G5)/N from Chapter 2 . Additionally, the critical 7 (3.3) causes the 


eigenvectors of H to be proportional to |s) ± |iu). As explained in Chapter 1, this 
causes the success probability to reach a value of 1. This is shown in figure |3.1 
for the cubic, cubic-quintic, and loglinear nonlinearities. A couple of observations 
are noteworthy. First, the cubic-quintic nonlinearity with one marked site has a 
wide peak in success probability, but with multiple marked sites, it has a narrow 
spike. Catching a narrow spike is more difficult than the wide peak, so searching 
with one marked site is “easier” than searching with multiple marked sites. This 
is counterintuitive, and it will be explicitly proven later. Second, for the loglin¬ 
ear nonlinearity, the success probability has a constant width. For the rest of the 


chapter, we choose 7 = 7 C as defined in (3.3). 
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f(p) = P, k= 1 


f(p) = p. k = 2 




f(p) = p - p 2 , k = 1 


f(p) = p-p 2 . k = 2 




f(p) = log(p), k = 1 


f(p) = log(p), k = 2 




Figure 3.1: Success probability as a function of time for search using the cubic, 
cubic-quintic, and loglinear Schrodinger equation with k = 1 and k = 2 marked 
sites and 7 at its critical value given by (3.3). The black solid curve is N — 100 
and the red dashed curve is N = 1000. The nonlinearity coefficient g scales as 
0(N ), O(IV), and 0(\/N/ log N) for the respective nonlinearities so as to make 
the runtime constant. 
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3.4 Runtime 


To derive the runtime of the algorithm, we follow the same procedure given 


in Chapter 2, generalized for (3.1). We begin by expliciting writing out (3.2), 


which yields two coupled, first-order ordinary differential equations: 

^ = i | ['yk + 1 + gfa\ ol + 7 \fkfN — k 


df3 

d t 


= z 


— ka + [7 (IV — k) + 


We define three real variables x(t), y(t ), and z(t) such that 


x = \a\ 


(3.4a) 

(3.4b) 

(3.5a) 


y + iz = a/3*. (3.5b) 

Note that x{t) is the success probability, which we want to find. To do this, we 


want to decouple (3.4a) and (3.4b) for a single differential equation in terms of 


x{t) alone and then solve it. We begin by differentiating (3.5a) by utilizing (3.4a): 

d|a | 2 


dx 

d£ 


1 = «TT- + !T a * = ^Vky/N^kz. 

d t d t d t 


Solving for z, 


da; 


z = 


2 7 V^\/iWl d t ' 

Note that the critical 7 depends on x: 


J 'k) J \ N-k 


(3.6) 


7c = V ' 1 + 9 


so we can use (3.6) to eliminate 0 in favor of x and da;/dt. We can also find 


an expression for eliminating dz/dt by differentiating this, but note that 7 = y c 
depends on time. Its derivative is 

1 


dru 

dt 


9_ 

N 


1 f> I - 1 - f’ 

k Ja N-k Jf3 


dx 
dt ’ 


where we’ve defined in analogy to f a and fp, 
d f(p) 


f = 

J OL 


dp 


r k k 


and fp = 


d f{p) 


dp 


\ f >\ 2 

L N-k N-k 
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Then the derivative of (J3.6[) is 
1 


dz 


■1 9 


d t 2\[k\jN - k \ N [k 


rf (r ) + 


1 fl-x 


N — k \N — k 


/da :\ 2 1 d 2 a; 
\ dt / 7 dt 2 


(3.7) 


So now we can eliminate dz/dt in favor of x, dx/dt, and d 2 x/dt 2 . 

Now let’s differentiate (3.5b) by utilizing (3.4a) and (|3.4b), which yields 


d d(a/3*) da d(3* 

di (9+!z) = ~6T = dU + “1T = ~ 2lkz+ ' 


+ 


7 Vky/N — k{ 1 — 2a;) j , 


where we’ve used 7 = y c to calculate the 27 k coefficients. Matching the real and 
imaginary parts, we get: 


^ = - 2 7 kz 

dt ’ 


d^ 


— = 27 ky + ryVkV N — k( 1 — 2x). 

kli 


In the first equation, we can eliminate 2 using (3.6), which yields 


dy 

d£ 


k da; 
N-kdt' 


This integrates to 


V = 


k 


-x + 


k 


k 


-(x-1), 


N-k V N-k V N-k 

where the constant of integration was found using ?/(0) = \fk(N — k)/N and 
x(0) = k/N. Using this to eliminate y in the second equation and simplifying, 

d£ 


dt 7 V N - k 


(N - 2 Nx + k). 


Eliminating dz/dt using (3.7) and simplifying, we get 

f' + 1 f {^ 2 

k 


d 2 x 1 g 
dt 2 ~ yiV 


N — k" 


) + 2'y 2 k (N — 2Nx + k ), 


which is entirely in terms of x and its derivatives. Plugging in for 7 = j c , 


d 2 a; 


N 


dt 2 1 + g(f a ~ // 3 ) N 


-f' + 1 f' 

k Ja N-k j0 


da; 

dt 


1 + g(fa ~ fp) 


N 


k (N — 2Nx + k). 


(3.8) 
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So we’ve decoupled (3.4a) and (3.4b), yielding a second-order ordinary differential 
equation for x. To solve it, let h(x) = (dx/dt) 2 so that dh/dx = 2d 2 x/dt 2 . Then 
we get a first-order ordinary differential equation for h(x): 

Id h Ng 


-f'+ 1 f 

2d t l + g(f a -fp)N [k Ja N-k Jf) 


h 


1 + g{f a ~ fg) 
N 


k (N — 2Nx + k). 


Solving this with the initial condition h(x = k/N) = 0, we get 

4 k(x - 1 )(k- Nx) [1 + g (f a - f/ 3 )} 2 


h(x) = 


N 2 


Taking the square root and noting that dx/dt = ±a Jh(x), 


dx 

d^ = ±1 


4 k(x - 1 )(fc - Nx) [1 + g (fa - fp)f 
N 2 


(3.9) 


We can solve this using separation of variables and integrating from t — 0 to t and 
x = k/N to x, which yields 


t = 


N 


-dx. 


(3.10) 


2Vk 7*0=6/K 1 + 9</o - M y (1 - x)(Nx - k) 

This integral depends on the form of f(p ). If it is analytically integrable, we 
get an expression for t(x), which we invert for x(t). For example, for the cubic 


nonlinearity, f(p) = p. Then f a — fp — (Nx — k)/(k(N — k )), and (3.10) can be 
integrated to yield 


t = -\ 


Nk 
k + g 


tan 


-1 


VNk \/1 — x 
y/k + g\/Nx — k 


7r 


(3.11) 


which can be solved for a success probability of 

N + (k + g) tan 2 
x(t) = - 


k+g j 
N 1 


N + ^r(k + g) tan 2 


This reaches a value of 1 at a runtime of 

1 7T y/N 


k+g_x 
N L 




\7. '7 2 


(3.12) 
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and the success probability is periodic with period 2i*. These results agree with 
Chapter 2. 

Returning to the general nonlinearity, if we are only interested in the run¬ 
time f* and not the entire evolution of the success probability, then we can instead 
integrate from x = k/N to 1: 


jl r“_i_,/ 1 dx 

2 Vk J„,k/N 1 + s(/o - Id) y (1 — i)(iVi - k) 


Evaluating this for the cubic nonlinearity yields (3.12), as expected. 


(3.13) 


3.5 Time-Measurement Precision 


As shown in FIGs. 2.4 and 3.1, the spike in success probability may be 
narrow. To quantify it, let’s find the width of the peak at height 1 — e. 


If we can explicitly integrate (3.10), then we can use the result to End the 


width in success probability. For example, the cubic nonlinearity yielded (3.11), 
which we use to find the time at which the success probability reaches a height of 
1 — e. Then the width of the peak at this height is 


At 


/ N -1 

VNkVe 

/ k + g 

_Vk + g\/N(l -e)-k_ 


We are interested in how this time-measurement precision scales with N, but the 
inverse tangent makes it difficult to see. Instead, Taylor’s theorem can be used to 
show that it suffices to keep the first term in the Taylor series: 


A*(°> 


2 Nk 
k + g 


k{N — k) 


+ 0(e 3 / 2 ). 


If we define G = g/[k(N — k)\, this agrees with our result from Chapter 2. 

For a general nonlinearity, we can hnd the leading-order formula for the 
time-measurement precision Aby Taylor expanding the success probability 
around x — 1, which is a maximum so the first derivative there is zero, and using 
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(3.8) for the second derivative: 


x(t) = x(t*) + x'(t*){t — £*) + -—j-—^-(t — f*) 2 + 


1 + 0 


1 + a\x=l ffi\x=l) 


N 


k(N - k)(t - 


This reaches a height of 1 — e at times 
£* ± 


t 


N 


1 + g{f a \x=i - fp\x=i)J k(N — k)' 
So, the leading-order width of the peak is 

Af<°> = 2N 


1 + g(fa \x=l - Mx= l) V k ( N ~ k ) 
For the cubic nonlinearity, f a \ x =i ~ fp\x=i — 1/A:, so we get 

At < 01 = 2 ' V 


(3.14) 


1 + g/k V Jb(JV - k) ’ ^ 3 ' 15 * 

which agrees with our previous result. 

To attain this level of time-measurement precision, say we use an atomic 
clock with lV c i oc k entangled ions. Then the time-measurement precision goes as 
0(l/lV c i oc k) [31f j30j. So the number of atomic clock ions we need is inversely pro¬ 
portional to the required time-measurement precision. This, plus the log IV qubits 
needed to encode the TV-dimensional Hilbert space, gives the “space” requirement 
of our algorithm. The product of “space” and time, which preserves the time-space 
tradeoff inherent in naive parallelization, gives the total resource requirement. 


Now that we have formulas for the runtime (3.13) and time-measurement 


precision (3.14) for a general nonlinearity of the form (3.1), let’s calculate them 


for the specific examples of the cubic-quintic and loglincar nonlinearities. But for 
comparison’s sake, let’s first review the results for the cubic nonlinearity. 


3.6 Cubic Nonlinearity 


The cubic nonlinear Schrodinger equation has the form (3.1) with f(p) = p. 

7TA /N 


From Eqs. (3.12) and (3.15), we found 


yjk + g 2 
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and 


A t<°> 


2 A 

l + g/k 


both of which agree with Chapter 2. If g 
1), then these become 


/ k(N-ky 

0(N K ) and k = 0(N X ) (with 0 < A < 


O (N ~ K / 2+l / 2 ) , k > A 

o (AT-V2+1/2) ; K < A 


= 


and 

O (iv- K+A / 2+1 / 2 ) , k > A 
O (iV _A / 2+1 / 2 ) , k<A 
To achieve this level of time-measurement precision, the number of ions in an 
atomic clock that utilizes entanglement must scale as the reciprocal of the precision 
pTTl f30 j. Including the log A qubits to encode the iV-dimensional Hilbert space, 
the total “space” requirement S scales as 


S 


O (AP'-V 2 - 1 / 2 ) , k > A/2 + 1/2 
O (log A), k< A/2 + 1/2 


Then the total resource requirement is 


o (a k / 2-a / 2 ) , 


ST 


< o (A _k / 2+1 / 2 log A) , 


O (A _a / 2+1 / 2 log A) , 


k > A/2 + 1/2 
k > X,k < A/2 + 1/2 
K < A 


This takes a minimum value of ST = A _a / 4+1//4 log A = (N/k) 1 ^ 4 log A when 
k = A/2 + 1/2, and it makes the width At*- 0 ' 1 constant. 

Of course, the cubic nonlinear Schrodinger equation, or Gross-Pitaevskii 
equation, is an effective nonlinear theory that only approximates the linear evo¬ 
lution of the multiparticle Schrodinger equation describing Bose-Einstein conden¬ 
sates. As worked out in Chapter 2 for the case of a single marked vertex, and 
generalized here to multiple marked vertices, since Grover’s algorithm is optimal 
[5] for (linear) quantum computation, the number of condensate atoms A 0 must 
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be included in the resource accounting such that the product of the space require¬ 
ments and the square of the time requirements is lower bounded by N. That is, 

since there are IVo oracles, each responding to a log IV bit query, 

' 

o (jy- A /2+i/2 + N~ k+1 N 0 log TV) , k > A/2 + 1/2 


ST 2 = 


< O(N~ K+1 N 0 logN ), 


O (N~ X+1 N 0 log N) , 


k> \,k< A /2 + 1/2 =n(N). 
K < A 


Then 



k > A 
k < X 


In the first region, this bound is maximized when k — 1, corresponding to the 
constant-runtime solution and beyond which it doesn’t make sense to increase k. 
In the second region, it is maximized when A = 1, he., the number of marked 
vertices scales with N. In both of these cases, the bound takes its strongest value: 


( N 

No = n (logIV 

As expressed in Chapter 2, to the best of our knowledge, this is the first bound on 
the number of condensate atoms needed for the Gross-Pitaevskii equation to be a 
good approximation of the linear, multiparticle dynamics. 



3.7 Cubic-Quintic Nonlinearity 


The cubic-quintic nonlinear Schrodinger equation has the form (3.1) with 
/ (p) = p — p 2 . Then 


fa~ fp = 


-N(N - 2 k)x 2 + k(N 2 - kN - 2 k)x - k 2 (N - k - 1) 
k 2 (N — k) 2 ' 


Plugging this into (3.13), the runtime is given by an integral of the form 
Nk 2 (N - k) 2 r* =1 1 




2 y/k Jx 0 =k/N ax2 + bx + c y (1 — x)(Nx — k) 


dx, 


where 


a = - 9 N(N - 2k) 














b = gk(N 2 -kN- 2k) 
c = —gk 2 (N -/,• - !)- k 2 (N - k) 2 . 
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This is analytically integrable, and the solution is 


+ _ it Nk 2 (N - k) 2 V2 2 o + 6 +-2 a-b+y/K 

2 2 Vk y/Tjy/A _^£ + y /£(k- N ) \Jt-y/A(k-N) 

where 

A = b 2 — 4 ac 
E = a + b + c 


£ = 2a/c + 2ciV + 6(/c + A). 

Let’s find the scaling of this runtime when g = 0(N K ) and k = 0(N X ) (with 
0 < A < 1) by finding the scaling of the individual terms and putting them 
together. We have 

a — O ( N K+2 ) 


Then 


c = 


b = 0 (N K+X+2 ) 

C) (A K+2A+1 ) , k> 1 
O (iV 2A+2 ) , /c < 1 


A = 

J 

0 (_/v 2k+2A+4 ) 

1 

hi > 

0 


1 

0 (W K+2A_M 

l ). 


hi < 

0 


0 

(A K+A+ 2) , 

K 

> 

1 


S = i 

0 

(a k+a+2 ) , 

K 

< 

1,A 

< K 


0 

< 

(a 2A + 2 ) , 

K 

< 

1,A 

> K 


f 

0 

(A k+A+ 3) , 

K 

> 

1 


£= < 

0 

(A k+A+ 3) , 

K 

< 

1,A 

< hi 


0 

(A 2A+ 3) , 

hi 

< 

1,A 

> K 
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We also have 

2 a ~\~ b ~\~ y/~A 

This is different, however, from 


O (iV K+A+2 ) , k > 0 
O (N k/2+x+2 ) k < 0 


—2 a — b + y/A 


( 


0 (N K+2 ), 

K > 

0 (N K+2 ), 

o < 

0 (N x+2 ) , 

o < 

0 ( j\r K / 2+A+2 ) 

K < 


1 

K < 1, A < K, 
n < 1, A > k 
0 


because when k > 0, the dominant term in y/A is b, which cancels with —b. The 
expression 

£ + y/~A[k — N) 


is a little tricky. The dominant terms of £ and y/~A(k — N ) cancel in certain cases. 
That is, when k > 1 or n < 1 and A < n, then £ = 2 ak + 2 cN + 6(/c + N ) is 
dominated by the term. When k > 0, A = b 2 — 4ac is dominated by the 6 2 
term, so yfK{k — N ) is dominated by — 6iV. So in these regions, the bN's cancel 
out, and we should ignore it when computing £ + y/~A{k — TV), thereby making 
£ = 0{2ak + 2ciV + bk) and y/~A{k — N) = Oibk — ^jf(k — N )). If we add them 
together, we get 

2 ac 

£ + a/A(A; - TV) = 2ak + 2cN + 2bk - — (k - N). 

Note that 2 ak + 2cN + 2 bk is dominated by —2gkN 2 + 2k 2 N 3 , and —^{k — N) 
is dominated by 2gkN 2 — 2kN 3 . Adding these, the 2gkN 2 factors cancel, leaving 
£ + y/~A{k — N) dominated by 2k 2 N 3 . So 


£ + y/~A{k — N) — O {N 2X+3 ) . 


It’s easy to see (he., we don’t have to worry about cancellations) that the scaling 
is also N 2X+3 for other values of k and A. Combining our results, 

2a + 6 + VA Jo (N K+1 / 2 ) , k > 0 
^£ + y/A(k - N) \o(AT«/ 2 +V 2 ) , K<0 
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The expression £ — y/A(k — N ) is different (easier) because the dominant term no 
longer cancels. So we have 


£ - VA (k - N) 


O (iV K+A+3 ) , K > 1 
O (jV K+A+3 ) , 0</c<l,A</c 
O (iV 2A+3 ) , 0</c<l,A>/c 

O (N 2X + 3 ) , k < 0 


Then 

O (iV K / 2 - A / 2+1 /2) ) K > i 

—2a — 6 + VA _ O (A^ k / 2_a / 2+1 / 2 ) , 0 <k<1,A<k 

y^£ - VA(& - TV) O (A^ 1 / 2 ), 0</c<l,A>/c 

O (iV K / 2+1 /2) ? K < 0 

Then 


2a + b + \/A —2a — b + VA ( O (N K+1 ^ 2 ) , k > 0 

y^£ + VA(/c — TV) VA(fc-iV) \o(A K / 2+1 / 2 ) , k<0 

We also have 


VsVa 


(9 (j\T3K/2+3A/2+3j ? K > 1 

O (jV 3K/2 +3A/2+3) ? 0</C<l,A</C 

O (iV K+2A+3 ) , 0 < « < 1, A > K 

O (j\T K/2 +2A+3^ ^ K < 0 


Putting all this together, 


O (iV K / 2+1 / 2 ), k > 1 
O (AT-K/2+ 1 / 2 ), 0</c<l,A<« 
O (AT- a / 2 +V 2 ) , 0<k<1,A>k 
O (iV- A / 2+1 / 2 ) , K < 0 
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Note that our formula can be reduced to two cases. When k > 1, then A < k since 
0 < A < 1. Similarly, when k < 0, then A > k. So we have 




O (AT-«/2+l/2) ) A < K 

0(N~ X / 2+1/2), A>k 


which is the same runtime order as search with the cubic nonlinearity. 

For the time-measurement precision, note that f a \x=i ~ fp\x=i — (k — l)/k 2 . 


Plugging this into (3.14), the width of the spike in success probability at height 
1 — e is 

2 N 


At<°) = 


-e. 


1+ g(k-l)/k 2 V k(N — k) 

When k — 1, the g term disappears. So varying g, while changing the runtime, 
doesn’t affect the width. When k ^ 1, it is 


At< 0 > = 


2 N 


1 + g/k I/ k(N - k ) 


e, 


which is the same as the cubic nonlinearity’s formula. Putting these together and 
letting g = 0(N K ) and k = 0(N X ) (with 0 < A < 1), we get 


A*(°> = 


O (A 1 / 2 ) , A = 0 

O (jV— k+a/ 2 +i/ 2 ) , A ^ 0, A < K 
o (AT— a /2+i/2) , A ^ 0, A > k 


So the runtime of search with the cubic-quintic nonlinearity scales identi¬ 
cally to search with the cubic nonlinearity. Furthermore, when there are multiple 
marked sites, the time-measurement precision also scales the same. But when there 
is a single marked site, the time-measurement precision scales as N 1//2 ; which is the 
same as Farhi and Gutmann’s linear algorithm |3j. This distinction between single 
and multiple marked sites is evident in figure 3.1| So for a single marked site, all 
the speedup that comes from the nonlinearity can be utilized without the expense 
of increasing the time-measurement precision. Thus search with the cubic-quintic 
nonlinearity is able to achieve a jointly-optimized runtime and time-measurement 
precision of 0(1). 
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As explained in Chapter 2, Grover’s algorithm is optimal [5], so there must 
be additional resources such that the product of the space requirements and the 
square of the time requirements is lower bounded by N. For the cubic-quintic 
nonlinearity, say the physical system is a periodic array of N waveguides, each 
long enough that the electromagnetic wave propagating through it performs the 
calculation. So the length of the waveguide would scale with the runtime £*. 
Keeping the cross sectional area of the waveguide constant, the number of atoms 
in a waveguide would also go as £*. Since we have N waveguides, the number of 
atoms would go as AT*. If the runtime is constant, then the number of atoms goes 
as Q(N), satisfying the optimality proof’s lower bound. 

The amount of energy, or number of photons, can also be included in the 
resource accounting. Say a waveguide needs P photons in the incident beam for 
it to behave like Kerr media with quintic corrections. Then we would need PN 
photons for the whole array. But it’s reasonable to say P is constant, so this would 
scale as N, again satisfying the optimality proof’s lower bound. 

We would also need charge to create an electric held at the marked wave¬ 
guides. Say this takes a constant number of resources. There are k marked wave¬ 
guides, so the resources for this would scale as k. While this may scale less than 
N, the other physical resources already satisfy the optimality proof’s lower bound. 

These resources may seem excessive, but if they scale linearly with N , then 
it is no different than any other database that requires the N items in the database 
to be physically written somewhere. 

For other physical systems that are effectively described by the cubic- 
quintic nonlinear Schrodinger equation, such as Bose-Einstein condensates with 
higher-order corrections [33], the additional resource to the runtime and time- 
measurement precision is some number of particles N 0l each of which responds 
to a log A" bit query. As proved above, when there are multiple marked vertices, 
the cubic-quintic nonlinearity solves the unstructured search problem in the same 
way as the cubic nonlinearity. Then the lower bound Nq = Q(N/ logN) from the 
cubic nonlinearity carries over. With a single marked vertex, the cubic-quintic 
nonlinearity requires a constant number of atoms in an atomic clock to achieve 
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Success Probability 


Figure 3.2: The integrand of the runtime integral (3.16) is the solid black curve, 
and the red dashed and green dotted curves are the integrands of the lower bound 
integrals (3.17), all for N = 1024, k — 5, and g — 1. 


the necessary time-measurement precision, and this yields the same bound. Thus 
the strongest lower bound on the number of particles is the same as for the cubic 
nonlinearity: 

N 0 = n 

To the best of our knowledge, this is the first bound on the number of particles 
needed for the cubic-quintic Schrodinger equation to be a good approximation of 
the linear, many-body Schrodinger equation. 



3.8 Loglinear Nonlinearity 


The logarithmic nonlinear Schrodinger equation has the form (3.1) with 
f(p ) = logp. Then 


fa~ ft3 = log 


N — k x 
k 1 — x 


Plugging this into (3.13), the runtime is given by the integral 


jl r _i_,/ 1 dx 

1-Jk J n .t/N 1 +9l°g( i V i Ti) V (! -x){Nx - k ) 


(3.16) 


Although it’s unclear how to directly integrate this, it is possible to bound it. Let’s 
begin with the lower bound. Splitting the region of integration into two parts, the 
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runtime is bounded below by 
N 

U > - 7 


r 1/2 

1 

_ Jk/N 1 + g log 

(N-k 1/2 \ 

l k 1-1/2 J 


-dx 


' 1/2 i + ^iog(^ T ^)y {i-x){n -i-k) 


dx 


These integrands are shown in figure 3.2 along with the original integranc 
trating that they are indeed lower bounds. These integrate to 


U > 


N 


2 Vk VN -k 


2(N-2k ) 


N l + gl 0g ^) 


,1 /.V - (^ + g\og(^^j 

"■ ■ / E x 1 v ' 


e 2 s - 


2.9 


9 V k 

where E] is the exponential integral 

r°° e -* 

Ei(x) = J —d t, 

which is bounded by 

^e _a log fl + < Ei(x) < e~ x log (l + - 

2 \ x J \ x 

Then the runtime is lower bounded by 
N 1 


t, > 


2 y/k VN -k 


2{N -2k) 


N i + ^i 0 g(hf_fc) 


log 1 + 


2 9 


\ 1 + 9 log ( 2(iV fc k) 

Now assume that g = 0(N K ) with k > 0. Then for large N, this becomes 

t* = n 


N 1 


k a log(f) 


(3.17) 


, illus- 


Now for the upper bound, we can again split the region of integration into 
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Success Probability 


Figure 3.3: The integrand of the runtime integral (3.16) is the solid black curve, 
and the red dashed and green dotted curves are the integrands of the upper bound 
integrals (3.18), all for N = 1024, k — 5, and g — 1. 


two parts: 


£* < 


N 

2 Vk 


" 1/2 


lk/N 


(1-1/2 )(Nx-k) 


dx 


(3.18) 


'i /2 1+jlog(tfc4_i£_ ^ (l-x)(N/2-k) 


dx 


These integrands are shown in figure 3.3 along with the original integrand, illus¬ 
trating that they are indeed upper bounds. The first region, however, is a poor 
bound, so we expect our result to not be tight. These integrate to 


t* < 


N 

2 Vk 


2\/N - 2k 
N 


+ 


VN^2k[l + glog(^)] 


Again assuming that g = 0(N K ) with k > 0 and large N, 


U = 0 I A I- 


But this does not provide much insight. It simply says that the nonlinear algorithm 
is no worse than the linear algorithm. This is expected because our upper bound 
is not very tight. 
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Success Probability 


Figure 3.4: Plot of (3.19) for N = 1024 and k = 5. The black solid curve is 
original logarithm, the red dashed curve is the bound from k/N < x < 1, and the 
green dotted curve is the bound from 1/2 < x < 1. 



Figure 3.5: The integrand of the runtime integral (3.16) is the solid black curve, 
and the red dashed and green dotted curves are the integrands of the upper bound 
integrals (|3.20|), all for N = 1024, k = 5, and g — 1. 
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To find a tighter bound for the runtime, we instead replace the logarithmic 


term in the denominator of the runtime integral (3.16) with a smaller function. In 


the region k/N < x < 1/2, we can use the line connecting those points, and in the 
region 1/2 < x < 1, we use the first-order Taylor approximation at x = 1/2: 


log 


N — k x 
k 1 — x 


< 


f= 2 ifc 1°S (nrO (Nx-k) k/N <x <1/2 

tog(^) +4(x-i) 1/2 < x < 1 


(3.19) 


The bounds for this logarithm are shown in figure |3~4| Then the runtime is bounded 
by 


t. < 


N 


,1/2 


2 Vk Jk/N i + S'/vzW lo s (nr) ( Nx ~ k n - x )( Nx ~ k ) 
N f 1 1 


dx 


2 y/k J 1/2 1 + 9 (log (^) + 4 (x - I)) V (1 - x){Nx - k ) 


dx. (3.20) 


These integrands are shown in figure 3.5 along with the original integrand, illus¬ 
trating that they are indeed upper bounds. They are also tighter than our previous 
attempt illustrated in figure |3.3[ The runtime integrates to 

-2 y/N - 2k 


t* < 


N 


l 


2\[k \ _2 k + 2g(N - k) log (S) 

Vn 


x tan 


+ 


7r 


y/N-2k + 2g(N-k) log(^) 

N 2 - 3kN + 2k 2 


yfW^kyfN V N-2k + 2g(N - k) log (^) 


N-k ' 
k , 


9 tan -1 \/ 4 s fc + jv '- 2 s- ;v +s jv lo s( 

V A/iV=2fc^l+2g+ 9 log(^) 


^1 + 2g + //log (^) + iV - 2gN + g IV log (^) j 

For g = 0(N K ) with k > 0 and large IV, this is dominated by the second term and 
becomes 

„( [N 1 
t* = o 



9 l ° g(f), 
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N 


Figure 3.6: The runtime of search with the loglinear nonlinearity for k = iV 1 / 4 
and g = AT/ 8 / log (N/k). The black circles were numerically calculated from (3.16) 
for N = 500 000 to N — 1 000 000 with intervals of 10 000, and the red solid line is 
the best-fit curve £* = 1.226 jV 0 261 . 


Combining this with the lower bound, the runtime is bounded between 

M 1 < t <M _I_ 

— v * vW(¥)’ 

where the notation f\(N) < / 2 (iV) denotes f\ (N) = 0(/ 2 (iV)), which implies that 
fi(N) > / 2 ( N) denotes fi(N ) = Q(f 2 (N)). Numerically, the actual runtime seems 
to be closer to the lower bound. For example, when k = 1V 1//4 and g = IV 1 / 8 / log N, 
the regression shown in figure [3dj| yields a runtime scaling of 0(1V°' 261 ), whereas the 
lower bound is O(iV 1 / 4 ) and the upper bound is 0(N 5 / 16 ) = 0(iV°' 3125 ). Given the 
frequent appearance of the ratio N/k, we define R = N/k. Also, the nonlinearity 
coefficient g appears with a factor of log A, so we now let g = 0(R a / log R) rather 
than 0(N K ) from before. Then the bounds are 

R 1/2 ~ a <R< R 1/2 ~ a/2 . 


For the time-measurement precision, note that f a \x=i ~ fp\x =l = so 


(3.14) says the width of the success probability is zero. But from figure 3.1, that 


can’t be right. This incorrect results arises because (3.14) was derived by Taylor 


expanding the success probability about its peak, but for the loglinear nonlinearity, 
the second derivative at the peak is negative infinity. 
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To get around this, we instead Taylor expand about a nearby point 1 — e. 


So we need the first and second derivative, which from (3.9) and (3.8) are 


da; 4k(1 — x)(Nx — k) [1 + g (f a — fp)\ 

Tt =± \ - W* - 


and 


d 2 x 


N 


dt 2 1 + g(f a ~ fp) N [k Ja N-k 


t/„ + T^-rf S 


da; V 


df J 


i + gift* - fp) 


N 


k (N — 2Nx + k) 


For the loglincar nonlinearity 


f a ~ fI -3 = log 


N — k x 
k 1 — x 


and 


k 1 


k Ja N — k J ^ x 1 — x a;(l — x) 

So near x = 1 — e for small e and large N, 


fa - fp\ x=1 _ £ ~ log ( ~ 


Still for large N, the first derivative is 


da; 

df 


x=l—e 


and the second derivative is 

d 2 a; 


df 2 


x=l—e 




kg 2 i°g 2 (fp 

N 


Then the Taylor expansion is 


x(t) « (1 — e) + 



(t - ti- e ) 


1 kg 2 log 2 (f i) 

2 N 


it - 


where is the time in which the success probability is 1 — e. Now let’s consider 
the time in which the success probability reaches a height of 1 — e/2, which is 
closer to the peak of 1. For small e, the first derivative of x{t) in this region 
is decreasing towards 0 because the success probability is approaching the peak 
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(where its derivative is zero). That is, for small e, we’re considering the region 
after the success probability’s inflection point. Then the width St = H-e /2 — t i_ e is 
a lower bound for the width At = i* — fi_ e / 2 , where t M is the time when the success 
probability is 1 (he., the runtime). Then the Taylor expansion becomes 

i e n \ , [ke (Nl\ 1 kg 2 log 2 (f ±) 

1 - 2 “ (1 - 0 + Viv 9loE m ~ 2-V- (,5i) ' 

This is a quadratic for Solving it and keeping the highest order terms, 


St 


_ [n 1 


^ 2 log 2 (fi) V k g\og (f-^) 

So the width of the success probability at height 1 — e is bounded by 

A t = Q 



k 5log (f 9 


Note that in figure 3.1, we chose g = 0(\/N/ log N) since k was constant, and 


it resulted in constant runtimes and widths. So this bound seems tight, and it is 
further evidence that the runtime £* is closer to its lower bound. 

To achieve this level of time-measurement precision in an atomic clock that 
utilizes entanglement, we need the number of clock ions to scale inversely with At. 
Also including the log A qubits to encode the TV-dimensional Hilbert space, the 
total “space” requirement S is 

S = O (R"- l/2 + log N) . 

Then the total resource requirement when cr > 1/2 is 

(1 + i? 1 / 2-0- log A) < ST < (iT /2 + R 1 / 2 ~ rT / 2 log a) . 

This is minimized when a = 1/2, yielding 

log A < ST < R}^ a log A. 


The upper bound equals the cubic nonlinearity’s total resource requirement. So 
the loglincar nonlinearity is at least as good as the cubic nonlinearity in reducing 
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the time-space resources. Given the numerical results from figure 3J3 the actual 
total resources seem closer to the lower bound. 

Of course, there must be additional resources such that the product of the 
space requirements and the square of the time requirements is lower bounded by N 
(5]. If the physical system (e.g., a Bose liquid) has N 0 particles, then each particle 
can be at any of the N vertices of the graph, which requires log N qubits for each 
particle. Then the “space” requirement S is No log N plus the number of clock 
ions to achieve the necessary time-measurement precision. That is, 


S — O (N 0 log N + i? 0 ' -1 / 2 ) 


for large N. Then 


R 1 ~ 2(T (N 0 logN + R a ~ 1/2 ) < ST 2 < R 1 ~ a (No log IV + i? CT ” 1/2 ). 


Since this must be lower bounded by N, 

R 1 ~ 2<t (Nq log N + R a ~ 1/2 ) = fl(N). 


When a < 1/2, this bound is satisfied regardless of Nq. When a > 1/2, then 


No = n 


NR 


'2(7—1 


log IV 


As er increases, this bound also increases. But there is no reason to increase a 
beyond 1/2, at which N 0 = Q(N log N), because that gives the optimal product 
of space and time when ignoring A/), and numerically gives constant runtime. So 
we’ve given a quantum information-theoretic bound for the number of particles 
needed for the logarithmic nonlinear Schrodinger equation to describe the phys¬ 
ical system (e.g., the number of atoms in a Bose liquid), and to the best of our 
knowledge, it is the first such result. 


3.9 Conclusion 

Our results indicate that a host of physically realistic nonlinear quantum 


systems of the form (3.1) can be used to perform continuous-time computation 
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faster than (linear) quantum computation. In particular, we’ve quantified this 
speedup by analyzing the quantum search problem, and the particular choice 
of nonlinearity gives rise to different runtimes, requires different levels of time- 
measurement precision, and necessitates a different number of particles for the 
nonlinearity to be an asymptotic description of the many-body quantum dynam¬ 
ics. 

Chapter 3, nearly in full, is a reprint of the material as it appears in “Quan¬ 
tum Search with General Nonlinearities” in Physical Review A 89, 012312 (2014). 
D. A. Meyer and T. G. Wong both contributed significantly to the work. 



Chapter 4 


Quantum Search on Strongly Regular Graphs 


4.1 Introduction 


As explained in Chapter 1, the quantum search problem can be formulated 
as a quantum random walk on the complete graph of N vertices, as shown in 


figure 1.1 , where the randomly walking quantum particle is initially in an equal 
superposition over all N vertices. Then the system evolves in a two-dimensional 
subspace spanned by the marked vertex and superposition of non-marked vertices. 
The next step in difficulty would be search on a graph where the system evolves 
in a three-dimensional subspace, namely that spanned by the marked vertex, the 
superposition of vertices adjacent to (or “one away” from) the marked vertex, and 
the superposition of vertices non-adjacent to (or “two away” from) the marked 
vertex. But this is precisely a strongly regular graph, examples of which are shown 
in figure 4T] with the vertices of each respective subspace colored red, blue, and 
white. A strongly regular graph with parameters (N, k, A, /i) is a graph with N 
vertices, each with k neighbors, where adjacent vertices have A common neighbors 
and non-adjacent vertices have /i common neighbors. 

As one might expect, only certain choices of parameters (. N , k, A, /i) give 
rise to strongly regular graphs. One constraint is that the parameters satisfy 


k(k — A — 1) — (N — k — l)/i, 


(4.1) 


which is proved by counting the number of pairs of adjacent blue and white vertices. 


On the left hand side of Eq. 4.1 the marked red vertex has k neighbors, so there are 
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Figure 4.1: Some examples of strongly regular graphs. From left to right: the 
Paley graph with parameters (5,2,0,1), the Paley graph with parameters (9,4,1,2), 
and the square graph L 2 (3) with parameters (9,4,1,2). Without loss of generality, a 
“marked” vertex is colored red, vertices adjacent to it are colored blue, and vertices 
not adjacent to it are colored white. 


k blue vertices. Each blue vertex has k neighbors, one of which is the red marked 
vertex, and A of which are other blue vertices. So it is adjacent to k — A — 1 white 
vertices. So the number of pairs of blue and white vertices that are adjacent to 


each other is k(k — A — 1). On the right hand side of Eq. 4.1 we count the number 


of pairs another way, beginning with the number of white vertices. There are N 
total vertices in the graph, one of which is red and k of which are blue. So there 
are N — k — 1 white vertices. Each of these white vertices has is adjacent to /i 
blue vertices. So there are (N — k — l)/i pairs of blue and white vertices that are 


connected to each other. Equating these expressions gives Eq. 4.1. Note this is a 
necessary, but not sufficient, condition for a strongly regular graph to exist. 

Equation 4.1| also implies that that k, the degree of the vertices, must be 
lower bounded by a /N. That is, 


k 2 > k(k -\-l) = (N-k- l)/i. 


Then 

k = f }(VN). (4.2) 

Additional constraints on the parameters (N, k , A, /i) divide strongly reg¬ 
ular graphs into two types [44, 45] : 


























59 


1. Type I graphs, also called conference graphs, satisfy 

(N- 1)(A - p) + 2k = 0, 
which means (N, k, A, p) can be parameterized by 

N — At + 1, k — 2t, A = t — 1, and p = t. 

This parameterization reveals that 

k = Q(N), (4.3) 


which will be useful later. Furthermore, Type I graphs exist if and only if N 
is the sum of two squares (one of the squares can be zero, so N = 9 = 3 2 + 0 2 
is acceptable). 


A large number of Type I graphs are Paley graphs, where N is congruent 
to 1 mod 4. The two smallest Paley graphs, where N = 5 and N = 9, are 


shown in figure 4.1 


2. Type II graphs satisfy 

(N-l)(p~ A) -2k AT , 
a 

This condition is not sufficient for the existence of a strongly regular graph. 
That is, just because a set of parameters (N, k, A, /j) satisfies this equation 
does not mean a strongly regular graph exists with such parameters. There 
are, however, some parameter families that do exist, three examples of which 
we now discuss. 


A. Square lattice graphs, an example of which is shown in figure |4.1[ can be 
pictured as a square lattice of t 2 vertices, where vertices are connected if and 
only if they are in the same row or column. They are denoted L 2 {t) according 
to the parameterization 


N = t 2 , k = 2{t — l), A = t — 2, and p — 2. 


B. Latin square graphs are similar to square lattice graphs, except each vertex 
is given a symbol that only appears once in each row and column. Then 
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Figure 4.2: The Latin square graph L 3 (3) with parameters (9,6,3,6). 


vertices with the same symbol are additionally connected. An example is 
given in figure |4~2] They are denoted Lz(t) according to the parameterization 


N = t 2 , k — 3(t — 1), A — t, and /i = 6, 


with t > 3. 

C. Triangular graphs, an example of which is shown in figure |4~3| are denoted 
T(t) and are parameterized by 

N — -t(t — 1), k — 2(t — 2), A = t — 2, and p — 4, 

with t > 4. They are constructed by labeling each vertex with a different 
unordered pair of different numbers in the set {1, 2,... ,t}. Two vertices are 
connected if they have a common number. 

Note that each of these examples of existing parameter families have 

k = Q(VN), (4.4) 


meaning they reach the lower bound given by Eq. 4.2 


It seems likely that almost all strongly regular graphs are asymmetric, 
meaning their automorphism groups are trivial. While this has not been proved in 
general, it has been proved for Latin square graphs, which were introduced above 
and shown in figure 4.2 [TO]. Thus, they lack global symmetry, which is intuitively 
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Figure 4.3: The triangular graph T(4) with parameters (6,4,2,4). 


believed to be necessary for fast quantum search [TJ. In this chapter, we show this 
intuition to be false, i.e., that a randomly walking quantum particle on strongly 
regular graphs optimally [5] solves the quantum search problem in 0(a//V) time 
for large N. 


4.2 Setup 


The vertices serve as a N -dimensional computational basis, which we label 
{10), 11),..., \N — 1)} and reduce to a three-dimensional subspace because there 
are three types of vertices: the red marked vertex, the k blue vertices that are 
adjacent to the red marked vertex, and the N — k — 1 white vertices that are not 
adjacent to the red marked vertex. Let’s call the basis states corresponding to 
these |u>), |a) (for adjacent), and | b) (since we called the other basis state |a)), 
respectively. That is, 


w) 


fl\ 

0 

w 


a) 


Vk ^ w 



l 

y/N-k-l 


i x ) 

(x,w)££ 


( °\ 

0 

V 1 / 
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The system begins in the equal superposition of all vertices, which we can write in 
the three-dimensional {|zt?), |a), | 6 )} basis: 


\s) = 


y/N 

1 

Vn 

i 

Vn 

i 

Vn 




\ X^W X'/'W / 

^|iu) + Vk\a) + VN — k — 1| b) 


\ 


1 

Vk 

^VN — k — 1 j 


The system evolves by Schrodinger’s equation with Hamiltonian 


H = —yL — |ta)(w|, 


where 7 is the amplitude per unit time of the randomly walking quantum particle 
transitioning from one vertex to another, and L is the graph Laplacian which 
effects a quantum random walk on the graph [I], More specifically, L = A — D, 
where A i:j = 1 if (i,j) G £ (and 0 otherwise) is the adjacency matrix indicating 
which vertices are connected to one another, and Du = deg(f) (and 0 otherwise) 
is the degree matrix indicating how many neighbors each vertex has. In the case 
of strongly regular graphs, each vertex has degree k, so the degree matrix is a 
multiple of the identity matrix: D = kl. This is simply a rescaling of energy, so 
we can drop it without observable effects. Then the Hamiltonian is 


H = —7 A — |tc) (w 


Now let’s determine H in the three-dimensional {|iu), |a), | 6 )} basis. The 
\w)(w\ term is simply a 3x3 matrix with a 1 in the top-left corner and 0’s everywhere 
else: 


w){w 


A 0 °\ 
0 0 0 
y0 0 Oy 
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Figure 4.4: The Petersen graph, which is a strongly regular graph with parameters 
(10,3,0,1). Without loss of generality, the vertices have been numbered and the 
marked vertex was chosen to be 0. 


The adjacency matrix A can be a little tricky, so let’s work it out explicitly for the 


Petersen graph, as shown in figure 4.4, in its (N = 10)-dimensional computational 


basis. With the labeling in figure 4.4, the adjacency matrix is 


A 


/o 1001100 
1 0 1 0 0 0 1 0 
0 1 0 1 0 0 0 1 
0 0 1 0 1 0 0 0 
1 0 0 1 0 0 0 0 
1 0 0 0 0 0 0 1 
0 1 0 0 0 0 0 0 
0 0 1 0 0 1 0 0 
0 0 0 1 0 1 1 0 
\0 0001011 


0 

0 

0 

1 

0 

1 

1 

0 

0 

0 


0 ^ 

0 

0 

0 

1 

0 

1 

1 

0 

0 / 
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The basis states {|iu), |a), \b)} in the computational basis are 



Then the adjacency matrix acting on each basis state is 
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So the adjacency matrix in the three-dimensional (|iu), |a), |6)} basis is 


A 


0 ^ 
y/3 0 

i 0 ^ 

\ U V3 



This is symmetric, as expected. 

Let’s examine where these terms come from. The normalization factors of 
| w), |a), and | b) are 1, l/\/3, and l/\/6, respectively. By the definition of the 
adjacency matrix, we can write it as a matrix of normalization conversions multi¬ 
plied component-by-component with another matrix that describes the connection 
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between vertices in \w), |a), and \b): 

1 A 7i\ (O 3 0\ 

A= VS 1 1 0 2 , 

W $ 1/ V° 1 2 / 

where Vindicates component-by-component multiplication (as in Matlab syntax). 
That is, the (i, j)-th entry of the first matrix is the normalization factor of the jth 
subspace basis vector divided by the ith ( e.g ., the entry at (2,3) is the normal¬ 
ization factor of | b) divided by the normalization factor of |a), or l/v^6 divided 
by 1 /x/3), and the (i,j)-th entry of the second matrix is the number of vertices 
in the jth subspace basis vector that are adjacent to a single vertex in the ith 
subspace basis vector (e.g., the entry at (2,3) is 2 because there are two white 
vertices connected to a single blue vertex). 

This example of the Petersen graph reveals what A should be for a general 
strongly regular graph with parameters (N, k, A,/i). First, recall that the normal¬ 
ization factors for |w), |a), and j b) are 1, 1 /Vk, and 1/y/N — k — 1, respectively. 
Then the adjacency matrix is 



Now let’s determine the second matrix. First, there are k blue vertices that go 
into the red vertex. So we have 



Next, there is one red vertex, A blue vertices, and (k — A — 1) white vertices that 
go into a blue vertex (for a total of k , the degree of a vertex). So we have 
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Finally, there are /j blue vertices and (k — /i) white vertices that go into a white 
vertex, so 


/ 


71 = 


1 

a Jk 


y/N-k-1 

\VN-k-l 1 


1 

y/k 

1 

~h 

y/k 


. 1 \ 
Viv-/c-l 

y/k 


/o k 0 \ 

1 A k — A — 1 
yO ii k — /i 


Multiplying these component-by-component, 

/ 0 


A = 


a Jk 
Vk A 

q y/ N—k— 1 ( 


0 




Vo 

The adjacency matrix must be symmetric. Let’s check the two non-obvious terms: 

a /k ? a JN - k- 1 

k(k — A — 1) = (iV — k — 1 )/r. 


But this is precisely Eq. |4.1[ so they are equal. Furthermore, it says that both of 
these are equal to yjjljk — A — 1. So we have 

/ 0 0 \ 


71 = 


a Jk A yj~jjk — A — 1 

^ 0 y fjlsjk — A — 1 k — /i j 


Thus the Hamiltonian H = —yA — |rc)(w| in the three-dimensional {|w), |a), \b)} 
basis is 

a fk 0 


( I 


H = - 7 


\ 


a Jk A y/j2yJk — A — 1 

y 0 yJjJi's/k — A — 1 k — /i j 


(4.5) 


4.3 Critical Gamma from Eigenstate Overlaps 


For the complete graph, the critical 7 caused the eigenstates of the Hamil¬ 
tonian to be proportional to |s) ± |w), which caused the system to evolve from |s) 
to \w) in time n\jN /2. We want the energy eigenstates to have the same form, 
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Figure 4.5: For the Latin square graph with parameters (900,87,30,6), overlaps of 
the eigenstates of the Hamiltonian in Eq. 4.5 with |u>) and |s) and the gap between 
the two lowest energy eigenvalues for various values of 7 . Note that |(s|' 02 )| J is 
near zero and not visible. 


but there are three of them |^ 0 )j |^i), and |^> 2 ). As shown in figure 4.5, this isn’t 
a problem—the projection of the highest energy eigenstate |- 0 2 ) onto |iu) and |s) is 
small, so it doesn’t contribute significantly to the evolution of the system. So only 
the ground and first excited states matter, and there is a point near the middle of 
the plot where they are roughly proportional to |s) ± |w). This corresponds to the 
critical value of 7 . 

So we expect to find a value of 7 that causes 


|s) ± |ta) = 


1 \ , , Vk , 4 y/N-k-1 

-j= ± 1 )\w + -A= a) + A-—- 

Vn J Vn Vn 


I b) = 


Vk 

Vn 

VN-k-1 

Vn 


to be an energy eigenvector. So let’s find the eigenvectors of H and choose 7 so 
that they have the desired form. The characteristic equation for the eigenvalues x 
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(since A is already used as a parameter of the strongly regular graph) of —H /7 is 


0 = det I —-i/ — y/ 


7 


X 


a fk 


Vk 

A-X 


y/Jly/k — A — 1 


0 y fjly/k — A — 1 k — fi — x 

= (7 ~ x) ~ X ^ k _ ^ _ * A - 1)] - VA; Vk(k - fi - x) 

= -X 3 + (- + k - A - fj\ x 2 


+ 


7 

— (A; + A — /i) — A(fc — /i) + /i(k — A — 1) + k 

7 


X 


H—A(fc — /i)- fi(k — A — 1) — k(k — n) 

7 7 

This is a cubic equation of the form ay 3 + fey 2 + cy + d = 0. Solving it seems really 
messy, as does finding the eigenvectors. So we’ll need another approach. 


4.4 Perturbation Theory in the {|w), |a), 16)} Basis 

Recall that the critical 7 for the complete graph can be found using degen¬ 
erate perturbation theory. Let’s try the same approach for strongly regular graphs. 
The leading order terms of the Hamiltonian in Eq. |4.5 depend on the relative sizes 
of the parameters, so we break the problem into two cases: when k scales as N, 


and when k scales less than N (but still bounded by VN as in Eq. 4.2). 


4.4.1 Case 1: k — Q(N) 

When k scales the same as N, the leading order terms in the Hamiltonian 


from Eq. 4.5 


are 


H (0) = - 


7 


0 0 \ 
A y/Ji'/k — A — 1 

^0 y^Jly/k — A — 1 k — /i j 


/l 

7 

0 
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Clearly, | w) is an eigenstate of this with eigenvalue —1. It’s straightforward to 
show that 

1 . 1 


| r) = 






X'/'W 


y/N^l 


Vk\a) + VN — k — l|6)j , 


which is approximately |s), is also an eigenvector of H (°), but with eigenvalue —y/c: 

0 0 \ / 0 \ 

A y/Jiy/k - A - 1 ^ 

^0 yjJly/k — A — 1 k — n y — k — lj 

( 0 


tf(°)|r) = - 7 


/l 

7 

0 


7 


\[ N ~^1 


~7 


~7 

VF^T 


-7 A; 


\ 

Aa/A; + y/Jiy/k — A — 1a/ TV — k — 1 
\VkyJ-jl\Jk — A — 1 + (A; — /i)\/iV" — fc — ly 

/ 0 \ 

A^A) + vF(A; - A - 1) 

— k — 1 + (A: — /^a/F — A; — ly 

/ 0 \ 


\fk{k — 1 ) 

\k\/N — k — 1 y 


/ 


^F^l 


\ 


0 

\fk 

JN — k — 1 y 


= -7 fc l r ), 


where we assumed that iV is sufficiently large such that k ~ k — 1. We want 
|in) and |r) to have the same eigenvalue so that the perturbation will cause the 
eigenstates of the full Hamiltonian to be a linear combination of |w) and |r). So 
we want — 7 A; = — 1 , or 

7 ci = (4.6) 

So with this value of 7 , the success probability should go to 1 for large N when k 
scales larger than a/F. This is verihed in hgure 4.6 for a Type I graph. When k = 


0 (vF), as in our examples of Type II graphs, we expect the success probability 


to not reach 1, and this is verihed in hgure 4.7 
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Figure 4.6: Evolution in each subspa ce fo r search on the Paley graph with pa¬ 
rameters (101,50,24,25) with 7 cl in Eq. 4.6. The black solid curve is Kwl^)! 2 , the 


red dashed curve is |(a|^)| , and the green dotted curve is 
the success probability nears 1. 


As expected, 



Figure 4.7: Evolution in each subspace for search on the Latin square graph with 
m — 30 (so parameters (900,87,30,60) with 7 cl in Eq. 4.6. The black solid curve 


is |(?n|'0)|b the red dashed curve is |(a|?/^)| , and the green dotted curve is |(6| , ?/’)|" 
As expected, the success probability doesn’t reach 1. 
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Finally, the third eigenvector of H (0 ^ is 


|e 3 ) = _ (y'N — k — T|a) - Vfc|6)) , 

which is necessarily orthogonal to both |w) and |r) and has eigenvalue — (A — n)/k\ 


|e 3 ) = -7 


/ - 

7 

0 


0 


0 

A y fjly/k — A — 1 

Vo y/Jiy/k — A — 1 k — ii J 

l 0 


\ / 
1 


y[W ~^T 


o N 
Viv - fc -1 
V -sfk ) 


-7 


a/FF) 


~7 

TF^T 


-7 

^F^T 




A y/N — k — 1 — \fkyfji\jk — A — 1 
\\[J l yJk — A — l-v/F — A; — 1 — \fk{k — n)J 

( ° 

Aa/A" — fc - 1 — - fc - 1 

y \fk{k — A — 1) — y/k(k — fi) j 

( ° 

(a-//)Fv i 

V (A — fi + 1)(—vF) J 

( o \ 


-7(A-/i)- 


a/F^F 
= —7(A — /i)|e 3 ). 


a/ At — A: — 1 

V -a/s 1 


Since |iy) and |r) are degenerate eigenvectors of H^°\ degenerate perturba¬ 
tion theory says that the eigenstates of the perturbed system are linear combina¬ 
tions of them: 

IF) = otw\w) + a r \r). 

The coefficients a w and a r can be found by solving the eigenvalue problem 




73 


where H wr = ( w\H + H^\r) and 


H™ = 


-7 


/ 0 yfk o\ 

yfk 0 0 

v 0 0 Oy 


Evaluating the matrix components, we get 

1 - 7 ^ 


-7^ 


vt/V^T | I aw \ _ £ ' 0,1 


k \7v^T 

Since 7 = 7 c i = 1/fc, this is 

-1 

-1 


-1 


-1 


QV 




v'iva 

—1 I \ 01 


= E 4 


a r 


Of 7 , 


a r 


Solving this eigenvalue problem, we get eigenvectors 
1 

—=. with eigenvalue E + = — 1 + 

72 V 1 


Vn^i 


1 i\ 1 

—= with eigenvalue E_ = — 1- , 

V2\l y/N^l 


Then the eigenstates of 7/ are 


7±) = -j=(T\w) + \ r)) 


E+ = -1 ± 


with eigenvalues 

L± ~ ^v^I' 

Note that the energy gap is A E = / W_ , Since | r) ~ |s), the system evolves from 
|s) to nearly \w) in time £* = n/AE = ny/N — 1/2 « TiyfN /2. 

4.4.2 Case 2: k = o(N ) 

When k scales less than A, the leading order terms in the Hamiltonian from 


Eq. 4.5 


are 


H® = - 


7 


A 0 0 \ 

7 

0 A 0 

yo 0 k — /ij 
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Figure 4.8: Evolution in each subspace for search on the Latin square graph with 
m = 30 (so parameters (900,87,30,60) with y ' 2 in Eq. 4.7. The black solid curve 


is |(tc|' 0 )| , the red dashed curve is |(a-l^)| , and the green dotted curve is |(5| , i/>}|" 
As expected, the success probability doesn’t reach 1 . 


Clearly, the eigenstates of this are \w), |a), and | b) with corresponding eigenvalues 


- 1 , —yA, and —y(fc — fi). Although |r) « |s) isn’t an eigenstate, from hgure 4.7 


the system is hardly in \a) for the Latin square graph, so it might be sufficient 
to use | b) instead of |r). Doing this, we want | w) and | b) to have degenerate 
eigenvalues, or —q(/c — fi) — — 1 , which implies 


7 c2 = 




(4.7) 


where the prime denotes that this is different from another (better) critical value 
of 7 for the second case that we derive later. 

Figure [478] shows the evolution of search on a Latin square graph. It reveals 
that y ' 2 is worse than 7 cl = 1/k, even though the latter was already suboptimal. 
This is expected, however, since more terms were dropped in the leading-order 
Hamiltonian used to derive y ' 2 than 7 cl . It is not apparent which terms to 
drop in the Hamiltonian in this (|tc), |a), | b)} basis. 


4.5 Perturbation Theory in the {|w}, |r), 1 63 ) } Basis 

In the last section, we used perturbation theory in the (|rc), |a), | 6 )} basis. 
When k = Q(N) it was clear which terms were best to drop in the leading-order 
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Hamiltonian H^°\ When k = o(N), however, too many terms were dropped. In 
order to drop fewer terms, we switch from the {|ru), |a), | 6 )} basis we’ve been using 


to the basis {|u>), |r), |e 3 )}. We can transform the Hamiltonian Eq. 4.5 to this basis 
by multiplying T l HT , where 

(l 


T = 


(h> l r > l e 3>) = 


0 

V° 


0 

y/k 

WI 

y/N-k-l 

v/Fn 


0 \ 

y/N—k—1 

vtv-l 
Vk t 
vFR/ 


Note that = T J = T, and T is a reflection in the (|a), | 6 ))-plane about 
the line through the origin in the direction of the vector (\/l — £, y/l + £), where 
^ = \Jk/{N — 1). Multiplying T~ 1 HT, the Hamiltonian in the new basis is 


/ 


H' = 


-7 


k(N—2) 
N -1 


y/N^l 

Vky/N-k -1 -Vky/N-k -1 
V v/vm n- 1 


Vky/N-k -1 \ 

—v/fcViV-fc-l 

N-l 

(A-M)(JV-l)+fc 
iV-1 / 


(4.8) 


As before we break the problem into two cases: when k scales as N, and when k 
scales less than N. 


4.5.1 Case 1: k — @(7V) 

Although we’ve already solved this case, we can replicate the result in the 


{|w;), |r), |€ 3 )} basis. The leading-order terms in the Hamiltonian in Eq. 4.8 

0 \ 


are 




-7 


0 


/l 

7 

0 k 


0 


^0 0 A — l-i J 

Clearly, the eigenvectors of this are |w), |r), and |e 3 ) with corresponding eigen¬ 
values — 1, — 7 k, and —q(A — //), which is consistent with what we found in the 
{|iu), |a), 16)} basis. We want | w) and |r) to be degenerate so that the perturbation 
causes the eigenstates to be linear combinations of |w) and |r). The value of 7 
that does this is 


7ci = 


k' 


which is the same result as Eq. 4.6 
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4.5.2 Case 2: fc = o(N ) 

When fc scales less than N, which includes Latin square graphs where k = 
©(a /N), the leading-order terms in the Hamiltonian in Eq. 

/ 1 0 Vk \ 


4.8 


are 


H<® = 


-7 


0 k 0 
/fc 0 \ — \Xj 

It’s clear that |r) is an eigenvector of H ^ with eigenvalue —qfc. Then the two 
other eigenvectors have the form 

M 

0 . 

W 

We want one of these to have the same eigenvalue —qfc so that H ^ is degenerate: 

fc A / Cl \ 

= —qfc 


H® 


0 

\c 3 J 


Solving this in the two-dimensional subspace: 

1 y/k 


7 


Cl 


a Vk A — n) \c 3 


0 

VV 

= k 


det I 7 k 


Vk 

Vk A — fi — k. 


Cl 

= 0 


-fc (A — /i — k) — fc = 0 


7 


7= * (1 ' 


k — A + /i 


1 fc — A + /i 

fcfc — A + /i — 1 

1 / 1 

= t 1 + 


= H 1 + 


= T 1 + 


fc — A — 1 + fj, 

1 

(N-k- l)£+/x 
fc 


(JV-l)M 
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This is the critical 7 when k = o(N ): 


7 c 2 — 7 + 


k (TV- 1)/U 
Finding the corresponding eigenvector is straightforward: 


-7 


/ I 

7 

0 

Vk \ 


|AA 


M 

0 

k 

0 


0 

= — 7 A; 

0 

y\/A: 

0 

A -n) 


W 


w 

A - k 

7 

0 

\[k 

\ 

M / 


0 

0 

0 


0 = 

0 


y 0 x — /i — kj 


W 


W 


Using the third line (or equivalently the first line), 

Vkc\ T (A — yU — k)c3 = 0 

k — AT /i 
Cl = -77-c 3 


\fk 


Normalization requires c\ T c§ = 1: 


(A; A T /i) 2 1 2 


Co T Co = 1 


c 3 = 1 + 


(fc-A + fiy 

aT" 


2 \ - 1/2 


So the eigenvector is 


|c) = 


/ Cl \ 

0 

w 


= 1 + 


{k - A T nf 


2 \ - 1/2 


0 

V 1 ) 


Note that the third eigenvector of is 


IT 


k 


(k - A T t ) 2 


- 1/2 


/^TL\ 

0 

V 1 ) 


(4.9) 
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with corresponding eigenvalue 

(A - /i ) 2 + k(l - A + /i) 


—1 k — A -I - jJj) 


The perturbation 


HW = 


-7 


to ^ 0 \ 

k n k 

Vn u Vn 


V 


n_«_ 

u vW 


0 / 


causes the eigenstates of 7/7 0 ) + H d) to be a linear combination of |r) and |c): 


| i>±) = oi r \r) + a c \c) 


To find a r and a c , we solve the eigenvalue problem 


TTipip I Iy* g 

■ Her H cc 


(X,n, 




= E ± 
a r ) \a r 


where H rc = + H^ ] \c), etc. These terms are straightforward to calculate. 

We get 

- 7 k fa,,\ _ fa,. 

k 


= E. 


± 


a r 


a r 


where A = (1 + (k — A + /i) 2 //c) 1,/2 is the normalization constant of |c), and we 


used Eq. 4.1| for the off-diagonal terms. For large N, this is 

—7 A; 


N.. . . 

k l l \ I 


fl- L 


= E -4 


-qfc 


a r 




a r 


Solving this, we get eigenvectors 


1 / 1 \ IN 

-j= I I with eigenvalue E + = — 7 k — 7 Ay — n 


1 

71 



with eigenvalue EL 


—7 A; + 7 AI 











79 


Then the eigenstates of H + H T) are 


l^) = ^(|r)T|c)) 


with eigenvalues 


E± = - 7 k ± 7 HW —/i. 


Now let’s find the success probability as a function of time. The initial 
state of the system is |s), which is close to |r). So the evolution of the system is 
approximately spanned by these two eigenstates: 

where H' = H + /7 (1 ' is the approximate Hamiltonian. That is, we’re ignoring 

/o 0 o\ 

# (2) = -7 0 0 0 • 

V° 0 &/ 

Then the state of the system approximately evolves in the subspace spanned by 

| 4 ±>: 

I i>(t)) ~ e~* E+t \ip+)(ip+\s) + s). 

Note that (^±( 5 ) = ^ ((r|s) =F (c|s)) ~ ^(1 T 0) = Then the success ampli¬ 
tude as a function of time is 

MW)} « 4 (e- £+, («.|V+) + e-‘ E -‘{w\^)) . 

Now let’s compute the inner products (w |4±)- To evaluate them, note that 


k — A T /1 — k — A — lT/iTl 

= (N — k — 1) -—|-/i + l 

= 4 V “ 1) + 1 




(4.10) 
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where we’ve used Eq. 4.1, k = o(N ), and assumed large N. Then (w\^)±) = 
=F^(u'|c) ~ ^f\A/j,N/k 3 / 2 . Plugging this in, 

Mm) *= 1^2 (-e- E +‘ + e-‘ E ->). 

Plugging in for the energy eigenvalues, 

(wMt)) « e~ ilkt -A^- + jrtAyfmv* 

2 k 3 / 2 \ 

The exponentials give us 2i sin(-). 

(«#(t)) ~ e ~ iyktA ^/2 isin 

Then the success probability is 


\{wm ))\ 2 « ( A ^p\ shr IryAJjfit ) . 


Using 7 « 1/k, 




AgX*? (aY 


k 3 / 2 


A: 3 / 2 ' 


Using Eq. 4.10, the normalization constant of |c) becomes 


A= 1 + 


(k — A + p) s 


2 \ - 1/2 


1 + 


(fiN) 2 Y 1/2 k 312 

~ [iN ’ 


k J \ ' k 3 

when k scales less than or equal to (/ilV) 2 / 3 , which is true for the known parameter 
families of Latin square graphs (which are proved asymmetric [45]). pseudo-Latin 
square graphs, negative Latin square graphs, square lattice graphs, negative Latin 
square graphs, square lattice graphs, triangular graphs, and point graphs of partial 
geometries ea- For these, the success probability is 

|(w|V’(/))| 2 ~ sin 2 / 


VnJ ’ 

which is 1 at time £* = n\/N /2, which is the same as on the complete graph. This 


is shown in figure 4.9 for a Latin square graph, and it outperforms figures 4.7 and 


4.8, as expected. We expect the success probability to approach 1 for large N, 


which we confirm in figure 4.10 
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Figure 4.9: Evolution in each subspace for search on the Latin square graph with 
m = 30 (so parameters (900,87,30,60) with 7 c2 in Eq. 4.9 


The black solid curve 
2 


is |(w| , 0)| , the red dashed curve is |(a|-0)| , and the green dotted curve is |(6|^)| 
As expected, the success probability nears 1. 



Figure 4.10: Evolution in each subspace for search on the Latin square graph 


with m = 50 (so parameters (2500,147,50,6) with q C 2 in Eq. 4.9 The black solid 


curve is |(u>| , 0)| , the red dashed curve is | (a,IV 7 )I , and the green dotted curve is 
As expected, the success probability nears 1. 
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Thus, we’ve shown that quantum search on known strongly regular graphs 
behaves like search on the complete graph for large N, reaching a success proba¬ 
bility of 1 at time 0(y/N). This requires choosing 7 = 7 c i = 1 /k when k = ©(TV) 
and 7 c 2 = 1 /k + l/[(iV — 1 )/r] when k = o(N). Since this includes strongly regular 
graphs that are asymmetric, it disproves the intuition that fast quantum search 
requires global symmetry. 

Chapter 4 is based on a paper, “Global Symmetry is Unnecessary for 
Fast Quantum Search,” published in Physical Review Letters 112, 210502 (2014). 
J. Janmark, D. A. Meyer and T. G. Wong all contributed significantly to the work. 



Chapter 5 


Nonlinear Quantum Search on Sufficiently 
Complete Graphs 


5.1 Introduction 


In Chapter 1, we showed that a randomly walking quantum particle can 
locate a marked vertex on the complete graph with probability 1 in time n\/N / 2 . 
It does this by evolving in a two-dimensional subspace with energy eigenstates 


iV’op) 


J_ / Vn 
v^V Vn± 1 


(|s>±k», 


(5.1) 


at the critical 7 , where |s) is the equal superposition of the vertices, and \w) is 
the marked vertex that we want to find. The corresponding energy eigenvalues of 
these eigenstates are 


E 0 ,i — —It 


1 

7F 


(5.2) 


so the system evolves from |s) to \w) in time tt/AE = tt\/N / 2 . 

In the previous chapter, we showed that a randomly walking quantum par¬ 
ticle can search for a marked vertex on a strongly regular graph with the same 
asymptotic behavior as on the complete graph. That is, although strongly regu¬ 
lar graphs are not complete (and even lack global symmetry), they are “complete 
enough” for the search to primarily evolve in its two lowest energy eigenstates, 
which take the form (|r) ± \w))/\/2, where |r) is the equal superposition of non- 


marked vertices, at the critical 7 and for large N. This is the same as (5.1) up 
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Figure 5.1: The 4-dimensional hypercube, which has 2 4 = 16 vertices. A single 
marked vertex is colored red. Vertices “one away” (adjacent) are colored blue, 
“two away” are yellow, “three away” are green, and “four away” is white. 


to terms of order 1 /\/N. We call such graphs that evolve according to (5.1) with 
error terms that tend to zero for large N sufficiently complete. 

Strongly regular graphs are not the only sufficiently complete graphs—the 
hypercube is as well. The hypercube is even “less complete” than strongly regular 
graphs—whereas search on strongly regular graphs evolve in a three-dimensional 
subspace, search on the n-dimensional hypercube evolves in a (n + l)-dimensional 
subspace, which grows with N = 2 n . An example of this for four dimensions 


is shown in figure |5.1[ where vertices that evolve identically are the same color. 
Nonetheless, the hypercube is complete enough for search on it to behave like 
search on the complete graph. That is, search on the n-dimensional hypercube 


also primarily evolves in its two lowest eigenstates, which take the form of (5.1), 
up to terms of order 1/n, at the critical 7 and large N mmv\- As an exam¬ 


ple, the evolution of the success probability for the 10 -dimensional hypercube is 
shown in figure 5.2 while N = 2 10 = 1024 is large enough for the runtime to be 
near 7ta/1024/2 « 50.265, it is not large enough for the second peak to be near 
37r-\/l024/2 « 150.80, or for the success probability to be near 1 . 

In Chapters 2 and 3, we showed that a fundamental nonlinearity provides 
a computational advantage in searching on the complete graph. In this chapter 
show that nonlinearities sometimes provide the same speedup when searching on 
sufficiently complete graphs, depending on the nonlinearity and the graph. 
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Figure 5.2: Success probability as a function of time for search on the 10- 
dimensional hypercube, which has 2 10 = 1024 vertices, at the critical 7 derived 
by Childs and Goldstone [1]. 


5.2 Linear Search 


We begin by introducing notation to describe linear search on sufficiently 
complete graphs. First we assume the system exactly evolves with eigenstates 


and eigenvalues of the form (5.1) and (5.2), respectively, showing that it folds the 
marked vertex with probability 1 in time n\/N /2. Then we introduce error in the 
eigenstates and show how it propagates to the success probability and runtime. 
Since search on strongly regular graphs evolves in a three-dimensional subspace, 
and search on the hypercube evolves in a (n + l)-dimensional subspace, search on 
a general sufficiently complete graph may evolve in an M-dimensional subspace. 
That is, the N vertices of the graph can be grouped together in M sets m; of size 
\rrii\, where all vertices in a set evolve identically. Then the equal superpositions 
of identically evolving vertices 

1 


I = 


VW ■ 




mi 


form an orthonormal basis (|mo), \ m i)i ■ ■ ■, | ? 7 ^m-i)} for the M-dimensional sub¬ 
space. Without loss of generality, we pick the marked site to be | m 0 ) = |w), so 
|mo| = 1. In the case of a strongly regular graph with parameters (N, k, A, //), the 
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basis states of the subspace are 


|™o> 

|mi) 

| m 2 ) 



i x > 

(x,w)g€ 


1 

y/N-k-1 


J2 i x >> 

(x,w)g£ 


which correspond to the marked vertex, vertices adjacent to the marked vertex, 
and vertices not adjacent to the marked vertex. For the n-dimensional hypercube, 
we first label each of the N = 2 n vertices with an ri-bit string \z± ... z n ). Without 
loss of generality, we choose the marked vertex to be the string of all zeros: | w) = 
|0 ... 0). Then the vertices “one away” are bit strings with a single one (he., with 
Hamming weight 1), the vertices “two away” are bit strings with two ones (i.e., 
with Hamming weight 2), and so forth. Taking the superposition of vertices equally 
far from the marked vertex, we get 



Then the set { \rrik ) : k = 0,1,..., n} is an orthonormal basis for the (n + 1)- 
dimensional subspace. 

In this (|m 0 ), \mi ),..., \rriM-i)} subspace, the state |^(t)) of the system 
can be written as a linear combination of the basis states: 


M—l 


\m) = 5Z 


i =0 


We assume that the system evolves in its two lowest energy eigenstates, having the 


form (5.1). Then the amplitudes are 


Ci(t) = (rriil'ipit)) = (mi|e lHt \s) = (mi|- 0 o )(V ; o|s)e lEot + (m^)(^s)e lElt . 


Then from the dehnition of sufficiently complete graphs in (5.1), 


//iv 1 rw-^ 1 \ i / vn vn ±i 
{Ms) V2 ^Vn±i\ Vn) V2^Vn± 1 VN 
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and 


<m,|*, I ) = 4' / ^ U ' m< 


± 6, 


iO , 


y/2 V v^v± 1 y Viv 

where 5j 0 — 1 when i = 0 and 0 otherwise is the Kronecker delta. Then the 
amplitudes are 


it 


Ci(t ) = 


Jn cos \Jn) + 1 sm ( Jn 




i = 0 
i ± 0 


Squaring them, the probabilities are 


h(t)\ 2 = 


N 


b cos 2 -fe + sin 


vdv 


Jn) 


^ U) 


i = 0 
i / 0 


(5.3) 


This reveals that the success probability |co| 2 reaches 1 at time iry/N /2. 

Now, to introduce error, a sufficiently complete graph has eigenstates of 


the form (5.1) and corresponding eigenenergies of the form (5.2), both up to terms 


of some order e, where e tends to zero for large N. For strongly regular graphs, 
e = 1 /x/iV, and for the n-dimensional hypercube, e = 1/n. Propagating these 
errors through the previous calculations, we get probabilities 


\c 0 (t) r = 

and 


N +0U) 


COS 


VN 


+ O (e) ) t + [1 + O (e)] sin" 


VN 


+ O (e) ) t, 


\ciMt )\ 2 = 


| rrij 
N 


+ O (e) 


cos 


y/N 


+ O (e) ) t, 


where the Vs are inside the trigonometric functions. So if the error in the eigen¬ 


states (5.1) and eigenenergies (5.2) tends to zero for large N, then the evolution 


also tends to the complete graph’s evolution, as expected. 


5.3 Nonlinear Search 

For the nonlinear algorithm, we subtract from Hq an additional nonlinear 
“self-potential” V(t) = g/(|^(r, f)| 2 ), where / is a real-valued function, so that the 





















system evolves according to the nonlinear Schrodinger equation: 


d 


= [■ H ° ~.3/(IV’(M)l 2 )]V’(M)- 


V(t) 


Then for positive g this speeds np the buildup of probability amplitude. In the 
computational basis, the self-potential is 

N—l 


V(t) = g^2f{\(i\i >)\ 2 ) |«)(i|. 


(5.4) 


i=0 


Even with this nonlinearity, the system still evolves in the M-dimensional subspace 
spanned by {|m 0 ), |mi),..., \rriM-i)}- In this M-dimensional subspace, the self¬ 


potential (5.4) has off-diagonal elements equal to zero. Its diagonal terms are 

N -1 

(mi\V(t)\mi) =9^ Him?) 



o 

II 

= g 

1 

\m,i\ 


= g 

1 

\mi\ 



j&nii 


mi 


= gf 




For ease of notation, let’s define 


fi = f 


mi 


Then in the M-dimensional subspace, the nonlinearity is 

(fo 0 ••• 0 \ 

V(t) = g 

0 


\0 0 J'm~ i J 
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In the M-dimensional subspace, the equation of motion is governed by 

M—l 

H = -7 L - |m 0 )(m 0 | ~ g ^ 

*=o 

M—l 

= ~1 L ~ (1 + gfo) |m 0 )(m 0 | - g fi "7)("7 

2 = 1 

M—l 

= - 7 ^ - (1 + gfo ~ gfi ) |m 0 )(m 0 | - 3 / 1 I - g^(fi~ fi) \mi)(mi\. 

i=2 

The term proportional to the identity matrix can be dropped since it is a rescaling 
of energy (or an overall phase), which is unobservable. We want to show that there 
exists a critical 7 that causes the nonlinear evolution to follow the same path as 
the linear evolution for large N, but with rescaled time r. That is, we want to 
show that 2 — /1 can be dropped in comparison to fo — fi when we use the linear 
evolution, but with t —> t. Assume for the moment that we can do this. Then we 
have 

H = -7 L - (1 + gfo - gfi ) |m 0 )(m 0 | 
for large N. Then the critical 7 for the nonlinear algorithm is 

7c = 7 l (1 + 9fo - gfi ), 

where 7 ^ is the linear algorithm’s critical 7 . At 7 c , we have 
H = (1 + gfo - gfi) (-7 lL - |m 0 )(m 0 |), 

which is the linear Hamiltonian at its critical 7 with a rescaled factor that depends 
on fo — f\ ■ For large N, the leading-order behavior of f 0 — f\ is 

h~h= /(A cos 2 (^) + shr (^)) - / (A cos 2 (^)) 

which is the same for all sufficiently complete graphs that can be sped up by the 
nonlinearity, including complete graphs, and so we expect the nonlinearity to speed 
them up the same way. 

Now let’s prove that such a critical 7 exists, that we can drop /*> 2 — fi 
compared to fo — fi- This depends on the graph and the form of /. Assuming as 
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Figure 5.3: Success probability for search with a cubic nonlinearity and g = 
N — 1 on strongly regular graphs at the critical 7 . The black solid curve is 
with parameters (N,k,\,g,) = (509,254,126,127), and the red dashed curve is 
(4001,2000,999,1000). With this choice of g, the runtime is constant for large 
N. 


before that the error in the eigenstates (5.1) and eigenenergies (5.2) is of order e, 
let’s first work out /o — fi and /*>2 — fi for f(p) = p (i.e., the cubic nonlinearity) 
so that fi = |ci| 2 /|mj|. We get 


fo-fi = 0 (e) cos 2 


Vn 


+ O (e) r + (l + 0 (e)) sin 2 


Vn 


+ O (e) r 


and 


fi> 2 — fi — O 


ImA 


cos 


+ O (e) r, 


\rnfj ~~~ 

where the r’s are inside the trigonometric functions. Clearly, / 0 — f± dominates 
fi >2 — fi at later time because of the sine piece, but what about at short time when 
cosine dominates? Recall that the set mo corresponds to the marked vertex, and 
the other sets m 8 yo correspond to identically evolving vertices. Then for strongly 
regular graphs, the m^o’s correspond to the k = 0(v/7V) vertices adjacent to the 
marked vertex and the N — k — 1 = 0(iV) vertices not adjacent to the marked 
vertex. Since both |mi| and Imfl scale as N to a positive power, /j> 2 — fi scales 
smaller than / 0 — / 1 , so we can drop it. Thus a critical 7 exists that causes 
search with a cubic nonlinearity on strongly regular graphs to behave like search 


on the complete graph, as shown in figure 5.3 This argument doesn’t work with 
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Figure 5.4: Success probability for search with a cubic nonlinearity and g — N — 1 
on the hypercube. The black solid curve is the 10-dimensional hypercube with 
N = 2 10 = 1024 vertices, and the red dashed curve is the 11-dimensional hypercube 
with IV = 2 11 = 2048 vertices. 


the hypercube, however, because there is an m* that corresponds to the vertex 


furthest from the marked vertex, which is the white vertex in figure 5.1, so it has 
size \rrii\ = 1. For this vertex, /j > 2 — /i scales the same than fo~ fi, so we can’t 


use this argument to drop it. Nonetheless, figure |5.4| indicates that a critical 7 
might still exist for the hypercube—our argument has not ruled it out. 

These arguments persist when we consider nonlinearities where f{jp) is a 
polynomial, which includes the cubic-quintic nonlinearity. If the polynomial is 
order g, so that f(p ) = 0 (p 9 ), then / 0 — fi would be 

1 


O (e p ) cos 2p 


y/N 


+ O (e) r 


plus terms with sine. We compare this to 


fi >2 — fi — O 


iml 


I rail 


cos 


2 V 


Vn 


+ O (e) r. 


Again we have the same sufficient condition for the existence of a critical 7 , that if 
the sets rrqyo grow with N, then we can drop f t >2 — fi compared to fo — fi- This 


is shown in figure 5.5 for strongly regular graphs. Note that the peak in success 
probability is not as wide as for the complete graph in figure |3.1| this is expected 
because the width broadens as the success probability approaches 1. That is, as it 
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Figure 5.5: Success probability for search with a cubic-quintic nonlinearity and 
g = N — 1 on strongly regular graphs at the critical 7 . The black solid curve 
is with parameters (N,k,\,n) = (509,254,126,127), and the red dashed curve is 
(4001,2000,999,1000). With this choice of g , the runtime is constant for large N. 


approaches 1, / 0 goes to zero while f\ remains nonzero. Then the rescaling of time 
slows down the evolution, causing a broad peak. For strongly regular graphs, the 
error with which the success probability approaches 1 decreases as N increases, 
so the peak is wider for large N. For the hypercube, figure |5.6| indicates that a 
critical 7 may exist, although the quintic term seems to have magnified the errors 


compared to the cubic case in figure [5i4| Since its success probability only reaches 
0 . 8 , we don’t expect the peak to broaden at all. 

For the loglinear nonlinearity f(p) = log p, we similarly find that whether 
/o — fi dominates f t — f\ depends on whether the sets of non-marked, identically 
evolving vertices have sizes that increase with N or not. Since log(a) — log( 6 ) = 
log(a/ 6 ), we have 


fo~fi = log 


= log 


[jf + 0 (e)] cos 2 + 0 (e)) r + [1 + 0 (e)] sin 2 (^= + 0 (e)) r 


Jr +0 (h7d 


cos^ 


Vn + °^) t 


[1 + O(JVe)] cos 2 + 0(e)) r + [N + O(iVe)] sin 2 (^ + 0(e) 


1 + 


COS- 


Vn + °^) t 
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Figure 5.6: Success probability for search with a cubic-quintic nonlinearity and 
g = N — 1 on the hypercube. The black solid curve is the 10-dimensional hypercube 
with N = 2 10 = 1024 vertices, and the red dashed curve is the 11-dimensional 
hypercube with N = 2 11 = 2048 vertices. 



We also have 


fi >2 - fi = log - 


jf +0 

y \ m i \ 2 


cos 2 

V^v+°( e )) 

r 

jf + ° 

J m d, 

’ 

cos 2 

Jn +°( e )] 

T 


= log 


l + O (i^r 

V l m 4 
1 + ofi^i 

V \mi\ 


At later time, /o — /i clearly dominates /,;>2 — /i because of the sine piece. At 
earlier time, it is dominated by the cosine piece, so it reduces to 

1 + O(Ne) 


log 


l + O 


Ne \ ’ 

l m il ) 


For strongly regular graphs, e = 1 /VN and \m\\ = £l{yf~N), so Ne/\rrii\ = 0(1). 
Then / 0 — fi for small time is dominated by log[0(Ae)]. Similarly, /)> 2 — f\ is 
dominated by \og[0(Ne/\rrii\)\, which is smaller, so we can drop it compared to fo~ 
/i, showing there exists a critical 7 . This is shown in figure 5.7| for strongly regular 
graphs. The second “peak” is strange because of numerical error; the derivative of 
logo: at x = 0 is nonzero, which makes the nonlinearity highly susceptible to noise, 


as shown in figure |5.8[ where the evolution begins to vary wildly shortly after the 
first peak. 
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Time 


Figure 5.7: Success probability for search with a loglinear nonlinearity and g = 
\/N/\ogN on strongly regular graphs at the critical 7 . The black solid curve 
is with parameters (N,k, A,p) = (509,254,126,127), and the red dashed curve is 
(4001,2000,999,1000). With this choice of g , the runtime is constant for large N. 


For the hypercube, we have the same issue as before; the set correspond¬ 
ing to the vertex furthest from the marked vertex has size 1 , so we can’t justify 


dropping /, ; > 2 — /1 compared to f 0 — fi. As shown in figure 5.9, our formulation 
doesn’t yield a critical 7 ; if it did, the peak of the success probability would stay 
near 0.8 with the nonlinearity. 

So whether a sufficiently complete graph can be sped up by the nonlin¬ 
ear Schrodinger equation depends on the graph and the nonlinearity. Nonetheless, 
we’ve shown that even with a degree of noncompleteness, some nonlinearities speed 
up search on certain sufficiently complete graphs in the same way as on the com¬ 
plete graph for large N. 

Chapter 5 is preliminary work for a paper to be published. D. A. Meyer 
and T. G. Wong both contributed significantly to the work. 
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Time 


Figure 5.8: Evolution of search with a loglinear nonlinearity and g = y/N / log N 
on a strongly regular graph with parameters (N, k, A, /i) = (509,254,126,127). The 
black solid curve is |co(t)| J , the red dashed curve is |ci(f)| 2 , and the green dotted 
curve is |c 2 (t)| 2 ; they correspond to the marked vertex, vertices adjacent to the 
marked vertex, and vertices not adjacent to the marked vertex, respectively. 



Time 


Figure 5.9: Success probability for search with a loglinear nonlinearity on the 
10-dimensional hypercube. The black solid curve is g = 0 (linear), the red dashed 
curve is g — 0.01, the green dotted curve is g — 0.05, and the blue solid curve that 
wildly varies is g = 0.5. 








Chapter 6 


Conclusion 


6.1 Summary 


To summarize our results, although quantum mechanics, which is governed 
by Schrodinger’s equation, is linear, there are nevertheless many-body quantum 
systems whose effective evolutions are governed by nonlinear Schrodinger equations 
of the form 

r)'ih _ _ _ 

( 6 . 1 ) 


^4? = [Ho-gnm ^ 


dt 

where / is a real-valued function. For example, Bose-Einstein condensates H21II31 
in] can be described by the Gross-Pitaevskii equation m 115] , which takes the 


form of (6.1) with f(p) = p. Nonlinear Kerr media with defocusing corrections 
can be described with f(p) = p—p 2 . Bose liquids [33] may be described 


with f{jp) = logp, which also retains the separability of noninteracting subsystems 


We showed that the nonlinear Schrodinger equation (6.1) can be used to 


solve the unstructured quantum search problem faster than standard quantum 
computation, using the cubic, cubic-quintic, and loglinear nonlinearities as spe¬ 
cific examples. In some instances, however, this causes the success probability to 
suddenly spike, which requires a certain number of atoms in an atomic clock to 
achieve the necessary time-measurement precision to catch the spike. This consid¬ 
eration of time-measurement precision as a physical resource is new. Even with it, 
we jointly optimized the runtime and time-measurement precision to outperform 


Grover’s algorithm, indicating that evolution by (6.1) is takes fewer resources than 
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evolution by Schrodinger’s equation, assuming the nonlinearity is fundamental. 

Of course, the nonlinearity is not fundamental, but arises as an effective 
description of the underlying linear dynamics. Since Grover’s algorithm is optimal 
[5], the speedup must be at the expense of increasing the number of particles in the 
physical system. Taking this into account, we arrived at bounds on the number 
of particles needed for the systems to be effectively described by the nonlinear 
equations. These are the first such bounds, and we novclly determined them by 
quantum information-theoretic means. 

The quantum search problem can be formulated as search on the complete 
graph, which evolves in a two-dimensional subspace. The next level of difficulty 
is search in a three-dimensional subspace, which strongly regular graphs support. 
Although strongly regular graphs are not complete, they are “complete enough” 
such that search on them behaves like search on the complete graph for large N , 
which we novclly used degenerate perturbation theory to show. Since this includes 
strongly regular graphs that are asymmetric, it disproves the intuition that global 
symmetry is needed for fast quantum search. 

The hypercube is even less complete than strongly regular graphs, yet is 
still “sufficiently complete” for search on it to behave like search on the complete 


graph for large N. The nonlinear Schrodinger equation (6.1) speeds up search on 
certain sufficiently complete graphs in the same way as it sped up search on the 
complete graph, so our results from nonlinear search on the complete graph carry 
over to these graphs. 

Thus we have proposed a new type of quantum computer that utilizes phys¬ 
ically realistic nonlinearities to compute in continuous time faster than (linear) 
quantum computation, even when some non-completeness is introduced into the 
underlying graph. 


6.2 Future Directions 

To physically implement our nonlinear search algorithms, the underlying 
search graph would need to be physically encoded in our three-dimensional world. 
Thus the complete graph, strongly regular graphs, and the hypercube are not viable 
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graphs to search on for large N. Square and cubic lattices, however, would be, and 
linear search on them has already been considered [T]. Preliminary results indicate 


that search by the nonlinear Schrodinger equation (6.1) also yields a speedup on 
arbitrary dimensional cubic lattices. But these graphs have sets of identically 
evolving vertices of constant size, so our analysis in Chapter 5 must be refined 
before we can analytically determine the precise computational advantage on cubic 
lattices. This would give a concrete runtime that experimentalists can seek to 
achieve, giving them a way to test the viability of our scheme for nonlinear, analog 
quantum computation. 
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